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ABSTRACT 

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) 
and 34 Palomar-Green (PC) quasars within z ~ 0.3 observed with the Spitzer Infrared Spectrograph 
(IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified 
using six independent methods that span a range in wavelength and give consistent results within ^ 
±10—15% on average. This agreement suggests that deeply buried AGN invisible to Spitzer IRS but 
bright in the far-infrared are not common in this sample. The average derived AGN contribution in 
ULIRGs is -35-40%, ranging from - 15 - 35% among "cool" (/25//60 < 0.2) optically classified HIL 
like and LINER ULIRGs to ^^50 and —75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. 
This number exceeds ~80% in PC QSOs. ULIRGs fall in one of three distinct AGN classes: (1) 
objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) 
systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, 
but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH 
equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions 
in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out 
the possibility that these objects represent a physically distinct type of ULIRGs. A morphological 
trend is seen along the sequence (1) — (2) — (3), in general agreement with the standard ULIRG — 
QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages 
of merger evolution. However, the scatter in this sequence, including the presence of a significant 
number of AGN-doniinated systems prior to coalesence and starburst-dominated but fully merged 
systems, implies that black hole accretion, in addition to depending on the merger phase, also has a 
strong chaotic/random component, as in local AGN. 

Subject headings: galaxies: active - galaxies: interactions - galaxies: quasar - galaxies: starburst - 
infrared: galaxies 
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1. INTRODUCTION 

More than twenty years ago, Sanders et al. (1988a, 
1988b) proposed the existence of an evolutionary connec- 
tion between ultraluminous infrared galaxies (ULIRGs; 
log[L(IR)/Z/0] > 12)^ and quasars. Their imaging and 
spectrophotonietric data on ten local ULIRGs were in- 
terpreted to imply that ULIRGs are dust-enshrouded 
quasars formed through the strong interaction or merger 
of two gas-rich spirals. This merging process had long 
been suspected to lead to the formation of an elliptical 
galaxy (e.g., Toomre & Toomrc 1972) and subsequent nu- 
merical simulations have lent support to this idea (e.g., 
Barnes 1989; Kormendy & Sanders 1992; Springel et al. 
2005; Bournaud et al. 2005; Naab et al. 2006). ULIRGs 
have been found since then to be an important popula- 
tion in the distant universe, a major contributor to the 
cosmic star formation at 2: > 1 — 2 (see, e.g., reviews by 
Sanders & Mirabel 1996; Blain et al. 2002; Lonsdale et 
al. 2006). The evidence is growing that the majority of 
the high-z ULIRGs are fed by continuous gas accretion, 
and only about one third arc gas-rich ( "wet" ) mergers 
(Daddi ct al. 2007; Shapiro et al. 2008; Genzel et al. 
2006, 2008; Forster Schreiber et al. 2006, 2009; Genel et 
al. 2008; Dekel & Birnboim 2008; Dekel et al. 2009). 
However, the fraction of ULIRGs involved in mergers 
appears to increase steeply at higher infrared luminos- 
ity (e.g., submillimeter-selcctcd galaxies; Tacconi et al. 
2006, 2008). 

Considerable effort has been devoted in the past decade 
to understand ULIRGs and quasars in the local universe, 

where galaxy merging and its relation to starbursts and 
AGN can be studied in greater detail than in the distant 
universe. Our group is conducting a comprehensive, mul- 
tiwavelength imaging and spectroscopic survey of local 
ULIRG and QSO mergers called QUEST Quasar and 
C/LIRG i;volution STudy. QUEST has already provided 
crucial new insights into merger morphology, kinematics, 
and evolution: we now know that ULIRGs are advanced 
mergers of gas-rich, disk galaxies sampling the Toomre 
merger sequence beyond the first peri-passage (Veilleux 
et al. 2002). The near-infrared (NIR) light distributions 
in many ULIRGs, particularly those with AGN-like op- 
tical and infrared characteristics, show prominent early- 
type morphology (i?^/^ law; Wright et al. 1990; Scoville 
et al. 2000; Veilleux et al. 2002, 2006). The hosts of 
ULIRGs lie close to the locations of intermediate-size (~ 
1 - 2 L*) spheroids in the photometric projection of the 
fundamental plane of ellipticals, although there is a ten- 
dency for the ULIRGs with small hosts to be brighter 
than normal spheroids. Excess emission from a merger- 
triggered burst of star formation in the ULIRG hosts may 
be at the origin of this difference. 

NIR stellar absorption spectroscopy with the VLT and 
Keck has also been carried out by our group to constrain 
the host dynamical mass for many of these ULIRGs. 
The analysis of these data (Dasyra et al. 2006a, 2006b) 

Woodlawn Drive, Honolulu, HI 96822. 
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1 Throughout the paper, L(IR) and L(FIR) refer to the 8 - 1000 
Hm and 40 — 120 nm luminosities as defined in Sanders & Mirabel 
(1996). 



built on the analyses of Genzel et al. (2001) and Tac- 
coni et al. (2002) and revealed that the majority of 
ULIRGs are triggered by almost equal-mass major merg- 
ers of 1.5:1 average ratio, in general agreement with 
Veilleux et al. (2002). In Dasyra et al. we also found 
that coalesced ULIRGs resemble intermediate mass el- 
lipticals/lenticulars with moderate rotation, in their ve- 
locity dispersion distribution, their location in the fun- 
damental plane and their distribution of the ratio of 
rotation/velocity dispersion [vrot sin(i)/fT]. These re- 
sults therefore suggest that ULIRGs form moderate mass 
(m* - 10" Mq), but not giant (5 - 10 x lO" M©) el- 
lipticals. Converting the host dispersion into black hole 
mass with the aid of the Mbh " relation (e.g., Geb- 
hardt et al. 2000, Ferrarese & Merritt 2000) yields black 
hole mass estimates ranging from 10'' ° Mq to 10^-^ Mq, 
with slightly larger values in coalesced ULIRGs than in 
binaries. BH masses derived from similar data on a dozen 
PC QSOs agree with those of coalesced ULIRGs (Dasyra 
et al. 2007). A recent analysis of HST/NICMOS data on 
several of these PG QSOs appears to support this con- 
clusion (Veilleux et al. 2009; see also Surace et al. 2001; 
Guyon et al. 2006). 

QUEST has also provided new quantitative informa- 
tion on the importance of gas flows in and out of ULIRGs. 
Direct evidence for powerful galaxy-scale winds has been 
found in most ULIRGs (e.g., Rupke et al. 2002, 2005a, 
2005b, 2005c; see also Martin 2005 and Veilleux et al. 
2005), and the metal underabundance and smaller yield 
measured in the cores of these objects (Rupke ct al. 2008) 
point to strong merger-induced gas inflows in the recent 
past as predicted by numerical simulations (e.g., Barnes 
& Hernquist 1996; Mihos & Hernquist 1996; lono et al. 
2004; Naab et al. 2006). 

The last two crucial issues addressed by QUEST are 
the nature of the energy production mechanism star- 
burst or AGN - in ULIRGs and QSOs and the im- 
portance of dust extinction along the merger sequence. 
Early optical and NIR spectroscopy has revealed trends 
of increasing AGN dominance among ULIRGs with the 
largest infrared luminosity, warmest 25-to-60 fim color, 
and latest merger phase (e.g., Veilleux et al. 2006 and ref- 
erences therein), but these results are potentially biased 
by dust obscuration. Mid-infrared (MIR) spectroscopy 
with the Infrared Space Observatory {ISO) (e.g., Gen- 
zel et al. 1998; Lutz et al. 1998b; Lutz et al. 1999; 
Rigopoulou et al. 1999; Tran et al. 2001) has provided 
crucial new information on the energy source in ULIRGs, 
less affected by the effects of dust, although the number 
of objects in the sample was limited by the relatively 
modest sensitivity of ISO. 

The most recent progress in this area of research oc- 
curred with the advent of the Spitzer Space Telescope. 
The Infrared Spectrograph (IRS) Guaranteed-Time and 
General Observation (GTO and GO) programs have pro- 
vided a wealth of new information on the physical prop- 
erties of local ULIRGs (e.g., Armus et al. 2004, 2007; Far- 
rah et al. 2007; Desai et al. 2007; Hao et al. 2007; Higdon 
et al. 2006; Imanishi et al. 2007; Lahuis et al. 2007; Spoon 
et al. 2007; CAo et al. 2008). In the present paper, we re- 
visit a carefully selected subset of these data and combine 
them with GO-1 IRS data acquired by our QUEST pro- 
gram, to calculate AGN contributions to the total bolo- 



3 



metric luminosities, quantify the issue of ULIRG evolu- 
tion along the merger sequence, and to determine where, 
if at all, quasars fit in this picture. 

Our earlier papers in this series have focussed exclu- 
sively on the QSOs. In Schweitzer et al. (2006; Pa- 
per I), we showed that starbursts are responsible for at 
least ~30%, but likely most, of the FIR luminosity of 
PG QSOs. We argued in Netzer et al. (2007; Paper II) 
that both strong- and weak-FIR emitting sources have 
the same, or very similar, intrinsic AGN spectral energy 
distributions (SEDs). In Schweitzer et al. (2008; Paper 
III), we found that emission from dust in the innermost 
part of the narrow-line region is needed in addition to 
the traditional obscuring torus in order to explain the sil- 
icate emission in these QSOs. The present paper reports 
the results from our analysis of the continuum, emission 
line, and absorption line properties of 74 ULIRGs and 
34 QSOs. In Section 2, we describe the sample. Next, 
we discuss the observational strategy of our program and 
the methods we used to obtain, reduce, and analyze the 
IRS spectra, including the archived data (Section 3, Sec- 
tion 4, and Section 5, respectively). The results are pre- 
sented in Section 6 and tested against the evolution sce- 
nario of Sanders et al. (1988a, 1988b) in Section 7. The 
main results and conclusions are summarized in Section 
8. Throughout this paper, we adopt Hq = 70 km s^^ 
Mpc-\ = 0.3, and = 0.7. 

2. SAMPLE 

The basic properties of the ULIRGs and quasars in 
our sample are listed individually in Table [TJ For a sum- 
mary of the properties of the ULIRGs by spectral types, 
infrared colors and luminosities, and morphology, see Ta- 
ble [TlJ The ULIRG component of our program focuses 
on the 1-Jy sample, a complete flux-limited sample of 
118 ULIRGs selected at 60 /im from a redshift survey 
of the IRAS faint source catalog (Kim & Sanders 1998). 
All 1-Jy ULIRGs have z < 0.3. Twenty-nine objects 
were observed under our own Cycle 1 medium-size pro- 
gram (#3187; PI Veilleux; Note that the 1-Jy ULIRG 
Mrk 1014 is also PG 0157+001). These objects were se- 
lected to be representative of the 1-Jy sample as a whole 
in terms of redshift, luminosity, and IRAS 25-to-60 /im 
colors. These data were supplemented by archival IRS 
spectra of 39 other galaxies from the 1-Jy sample, and 5 
archival IRS spectra of infrared-luminous galaxies from 
the Revised Bright Galaxy Sample (RBGS, Sanders et 
al. 2003; these objects are UGC 05101, F10565-I-2448, 
F15250-h3609, NGC 6240, and F17208-0014). Most of 
the archival spectra are from GTO program #105 (PI 
Houck), and three are from GO program #20375 (PI 
Armus). These spectra cover bright sources in the 1 Jy 
sample, while ours are deeper exposures of fainter ones. 
Together, they represent almost 2/3 of the 1 Jy sample. 
The 5 RBGS spectra represent well-studied benchmarks 
from the local universe. Figure [T] shows the distribu- 
tions of redshifts, infrared luminosities, and 25-to-60 /im 
IRAS colors for the combined set of ULIRGs compared 
with that of the entire 1-Jy sample. We confirm that 
the ULIRGs in our study are representative of the range 
of properties of the 1-Jy sample. Optical spectral types, 
which are referred to extensively in this paper, are taken 
from Veilleux et al. 1999a and Rupke et al. 2005a for 
the 1 Jy sample, and Veilleux et al. 1995 for the 5 RBGS 



objects. 

The original QUEST sample of quasars has already 
been discussed in detail in Papers I and II and this 
discussion will not be repeated here. Suffice it to say 
that the original QUEST sample contains 25 z < 0.3 
quasars, including 24 Palomar-Green (PG) quasars from 
the Bright Quasar Sample (Schmidt & Green 1983) and 
another one (B2 2201+31A = 4C 31.63) with a B mag- 
nitude that actually satisfies the PG QSO complete- 
ness criterion of Schmidt & Green (1983). Nine other 
PG QSOs (PG 0050+124 = I Zw 1, PG 0804+761, 
PG 1119+120 = Mrk 734, PG 1211+143 = Mrk 841, 
PG 1244+026 [NLSl], PG 1351+640, PG 1448+273 
[NLSl], and PG 1501+106) observed under different 
Spitzer programs were later added to the quasar sam- 
ple. Figure O emphasizes the fact that the quasars in our 
study cover the low redshift and low B-band luminosity 
ends of the PG QSO sample, while Figure [3] shows that 
the ULIRGs and quasars in our study are well matched in 
redshift. Finally, note that two ULIRGs, Mrk 1014 and 
3C 273, are also PG QSOs; we treat them as ULIRGs for 
the purposes of this study. 

High-quality optical and NIR images obtained from the 
ground and with HST are available for all ULIRGs and 
quasars in the present sample (e.g, Surace & Sanders 
1999; Scoville et al. 2000; Surace et al. 1998, 2001; Guyon 
et al. 2006; Veilleux et al. 2002, 2006, 2009). In addition, 
high-quality optical spectra exist for all 1-Jy ULIRGs 
(e.g., Veilleux et al. 1999a; Farrah et al. 2005; Rupke 
et al. 2005a, 2005b, 2005c) and PG QSOs (Boroson & 
Green 1992), and a large subset of these objects also 
have been the targets of NIR JHK-band spectroscopy 
by our group over the years (e.g., Veilleux et al. 1997, 
1999b; Dasyra et al. 2006a, 2006b, 2007) as well as some 
L-band spectroscopy (e.g., Imanishi et al. 2006a; Risaliti 
et al. 2006; Imanishi et al. 2008; Sani et al. 2008). These 
ancillary data will be used for our interpretation of the 
Spitzer data in Section 6 and Section 7. 

3. OBSERVATIONS 

Galaxies from our own program (#3187; PI Veilleux) 
were observed in the IRS modules SL, SH, and LH, using 
staring mode (Houck et al. 2004). Together, these mod- 
ules cover observed wavelengths of 5 — 35 fj,m. The high 
resolution data at observed wavelengths of 10 — 35 /xm 
(SH and LH modules, with resolution R ^ 600) allow 
sensitive measurements of important atomic and molec- 
ular emission lines. 

For targeting, moderate-accuracy IRS blue peak-ups 
were performed on the targets themselves rather than 
offsetting from 2MASS stars. This peak-up method is 
justified given the compact MIR continua in these sys- 
tems (e.g., Soifer et al. 2000; Surace et al. 2006). 

For the four binary ULIRGs with nuclear sep- 
arations exceeding 3" (F01166-0844, F10190+1322, 
F13454-2956, and F21208-0519), a unique observation 
was made of each nucleus. However, in 3 of these cases 
(F10190+1322 being the exception), aperture effects due 
to the larger slit sizes of the long-wavelength modules 
allowed accurate measurements of only one of the two 
nuclei. 

The observational setup used for the archival IRS spec- 
tra is described in detail in Armus et al. (2007) and ref- 
erences therein. It is essentially the same as the one we 
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used for our own program so direct comparison between 
the two data sets is permissible. 

Some objects in our sample have full low-resolution 
spectra {i.e., including both the SL and LL modules, 
covering 5 — 35 /im). We have used only the high- 
resolution data for spectral line measurements (except 
for the [Ne VI] line, which falls in the SL module). The 
LL data was used primarily for checking flux calibration. 
However, when only SL-I-LL data was available, or when 
the high-resolution data was of low S/N, the full low- 
resolution spectrum was used in the continuum fitting. 

Spitzer proposal ID numbers and exposure times for 
each IRS module are listed for all galaxies in Tabled] 

4. DATA REDUCTION 

For the majority of QUEST sources, we started with 
BCD data processed by version S12.0 of the IRS pipeline. 
For the non-QUEST 1 Jy ULIRGs that were reduced 
at a later date, data from pipelines S12, S13, or S15 
were used. Comparisons among these pipelines show only 
minor differences that do not impact our measurements. 

The data were first corrected for rogue pixels using an 
automatic search-and-interpolate algorithm (which was 
supplemented by visual examination). For the SL mod- 
ule, background light was then subtracted by differencing 
the two nod positions. For the SL and LL modules, the 
data was extracted prior to coadding. For SH and LH, 
we coadded exposures for a given nod position prior to 
extraction. 

The one-dimensional spectra were extracted using the 
Spectroscopic Modeling Analysis and Reduction Tool 
(SMART; Higdon et al. 2004). The extraction apertures 
were tapered with wavelength to match the point-spread 
function for SL and LL data and encompassed the entire 
slit for SH and LH data. The correction/extraction pro- 
cess was iterated until we were assured that the majority 
of hot pixels had been removed. 

For SL, we combined the two one-dimensional nod 
spectra for the orders SLl and SL2 separately and then 
stitched the orders by trimming a few pixels from one 
or the other order. (We discarded SL3 because of flux 
discrepancies.) For SH and LH, we combined the two 
nods and then stitched the orders together by applying 
multiplicative offsets for each order that were linear in 
flux density vs. wavelength (effectively removing a tilt 
artifact from certain orders where necessary). 

We subtracted zodiacal light from the high-resolution 
data using a blackbody fit to the Spitzer Planning Ob- 
servations Tool (SPOT) zodiacal estimates at 10, 20, and 
35/zm. 

Finally, we matched modules in fiux to form complete 
5 — 35 fj,m spectra. Because ULIRGs are compact mid-IR 
sources, different modules in general agree well in fiux 
at wavelengths where they overlap. Where there was 
disagreement, we used available low-resolution spectra 
(SL -|- LL) to improve the zodiacal subtraction, since 
these spectra were sky-subtracted using simultaneous sky 
observations. Where this was not possible, we used small 
additive offsets. 

To check the flux calibration, we computed synthetic 
IRAS flux densities at 12 and 25 /xm using the Spitzer 
data by averaging the flux densities over the IRAS band- 
passes. At 25 fiui, the agreement with IRAS is excel- 
lent (Figure [4]). The median IRAS-to- Spitzer Rux den- 



sity ratio is 1.04, with a standard deviation of 0.3. At 
12 /im, most of the IRAS fluxes are upper limits, but 
the agreement is still decent {fl2^^ / fi-z^*^'^^ ~ 1-25 on 
average for sources with /12 > 0.1 Jy). The cause of 
the small discrepancy at 12/im is unclear, but may result 
from Eddington-Malmquist bias. 

In this paper, we adopt the Spitzer- derived 12 and 
25 ijm fluxes to avoid the use of IRAS upper limits, a 
particularly severe problem for the LINER and Hll-like 
ULIRGs of our sample. 

5. DATA ANALYSIS 

5.1. Emission Lines 

In Tables [3] and IH we list atomic and molecular 
emission-line fluxes measured from our spectra. Mea- 
surements were made with the IDEA tool in SMART; we 
fltted Gaussian profiles atop a linear continuum. Upper 
limits were determined by assuming an unresolved line. 
The available resolution allowed us to decompose close 
line blends, including the important [Ne V] 14.32 /im / 
[CI II] 14.36 /im and [O IV] 25.89 /im / [Fc II] 25.99 /xm 
lalends. 

5.2. Continuum and Dust Features 

A vitally important task was to properly model the 
sum of the blackbody continuum emission, which is punc- 
tuated by deep extinction and absorption features, and 
the full-featured small dust grain continuum. The pri- 
mary goal of this modelling was to accurately extract 
the fiuxes of absorbed PAH features, but we also gained 
useful information about the continuum. Because we are 
not concerned with detailed physics, we have chosen a 
simple, but robust and empirically-motivated, method. 

Before fitting, we measured narrow-band fiux densities 
(of width 3.3% of the central wavelength) at regularly- 
spaced intervals across the continuum, avoiding deep ab- 
sorption features. These are listed in Tabled 

To fit the MIR spectra, we used the IDL package de- 
veloped to model the blackbody and silicate emission of 
QUEST quasars (Paper HI). We refer to this paper for 
the basic fitting details. Here we describe some unique 
features necessary for fitting the spectra of PAH-strong 
and sometimes deeply-absorbed ULIRGs. Typical fits 
are shown in Figure [5] 

1. A modifled version of the Chiar & Tielens (2006) 
Galactic Center extinction curve was used. To 
the basic extinction proflle, we added extinction 
from the water ice plus hydrocarbon feature at 
5.7 — 7.8 /im. The proflle is taken from observa- 
tions of F00183-7111 (Spoon et al. 2004). In 
place of the broad silicate features at 9.7 and 18 
/im we substituted features that were empirically 
derived from the deeply absorbed and almost com- 
pletely PAH-free ULIRG, F08572-h3915. The sil- 
icate proflle of this galaxy provides a universally 
good flt to the ULIRGs in our sample. However, 
using the original Chiar & Tielens silicate profile, or 
one from a less deeply-absorbed system, yields poor 
fits at high optical depths in a number of deeply ab- 
sorbed systems. The strengths of the water ice -|- 
hydrocarbon feature and silicate absorption -I- over- 
all extinction curve were allowed to vary indepen- 
dently. Foreground and mixed dust screens were 
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both tried, and we found that foreground screens 
fit better. 

2. We fit three blackbodies to almost all spectra. 
These blackbodies represent a convenient param- 
eterization of the MIR continuum. Due to the ab- 
sence of wavelengths probing the hottest and cold- 
est dust, the temperatures of these blackbodies do 
not represent actual dust components. Nonethe- 
less, their values provide useful guidance in under- 
standing the shape of the continuum. For most 
sources, three is the minimum number of black- 
bodies that produce an acceptable fit. However, 
the use of two or four blackbodies significantly im- 
proved the fit in six cases. 

3. The temperatures, water ice -I- hydrocarbon ab- 
sorption, and extinctions of the blackbodies were 
in general allowed to freely vary (see results in Ta- 
ble [6|). However, in some cases the water ice plus 
hydrocarbon absorption and/or the extinction had 
to be fixed to zero in a particular component, when 
it was apparent that the fitted value was unphysi- 
cal. In a handful of cases, the temperature of the 
hottest component was also poorly constrained by 
the fit. Experimentation and by-eye examination 
suggests that these unconstrained temperatures are 
not much larger than 1000 K, to which we fixed 
them. The actual temperatures are probably in 
the range - 700 - 2000 K. 

4. PAH emission was modeled using the average MIR 
spectra derived from the SINGS program galaxies 
(Smith et al. 2007). Pure PAH templates were 
created by running the PAHFIT program (Smith 
et al. 2007) to extract the PAH emission features 
from the four average SINGS spectra (including the 
17 /im emission band; van Kerckhoven et al. 2000; 
Peeters et al. 2004). The PAH model underlying 
PAHFIT is a series of Drude profiles; the 6.2 and 
7.7 /xm features, which we discuss in §6.3, consist 
of 1 and 3 Drude components, respectively. Note 
that the templates output by PAHFIT consist only 
of PAH emission, without an associated continuum 
from star formation regions; this continuum is fit 
by the blackbodies mentioned above. We chose two 
of these templates (from Smith's 3rd and 4th aver- 
age SINGS spectra) that provided the largest range 
in the ratio of the 7.7 /im and 11.2 PAH fea- 
tures. The overall strengths of these two pure PAH 
templates were allowed to vary in the fit. 

The fit results, listed in Table [3 were not sensi- 
tive to the choice of pure PAH template. Only a 
small range 0.13 dex) in PAH 6.2/7.7 /im and 
7.7/11.2 iim ratios was allowed by the two PAH 
templates and most galaxies are dominated by one 
or the other template. We found that using tem- 
plates produced from SINGS spectra of galaxies 
with very low values of the 7.7/11.2 /im ratio (J. D. 
Smith 2007, pvt. comm.) yielded poor fits. Thus, 
the strong suppression of the 6.2 or 7.7 /im com- 
plex observed in galaxies that host low-luminosity 
AGN (Sturm et al. 2006; Smith et al. 2007) and in 
some star-forming regions (Hony et al. 2001) does 
not occur in ULIRGs. 



We found that allowing PAH extinction below the 
level of Ay ^ 10 did not significantly affect the 
fits (i.e., the results were basically indistinguish- 
able from the Ay = case). We thus left the 
PAHs unextincted in most cases. Adding larger 
amounts of extinction had the effect of raising 
the 7.7/11.2 /im ratio and almost always made 
the fits worse. However, in four cases we allowed 
the PAH features to be extincted at or above the 
Ay ~ 10 level because the fit was significantly 
improved. These four cases are as follows, with 
PAH extinction in parentheses (two values indicate 
two fitted components): F00397-1312 {Ay = 44), 
F01494-1845 {Ay = 10/14), F20414-1651 {Ay = 
0/22), F21208-0519:N {Ay = 11). These values 
correspond well with the effective continuum ex- 
tinction (Section 6.2) in these sources: Ay ^ 50, 
18, 18, and 15, respectively. We discuss PAH ex- 
tinction further in Section 6.3. 

5. For four galaxies (Mrk 1014, F07598-h6508, 3C 273, 
and F21219— 1757) we included silicate emission 
components, as described in Paper III. 

6. Due to the absence of pipeline error spectra, the 
mixing of data over a significant flux range, and 
the different dispersion among different IRS mod- 
ules, careful weighting had to be performed during 
the fits (see Paper HI for details). We tried three 
methods: (1) weighting based on the actual fluxes; 
(2) weighting based on a power-law fit to the fluxes; 
and (3) an average of the two. In almost all cases, 
the third method produced the best fits. However, 
in isolated cases we used one of the other weighting 
schemes if it was clearly superior. 

Along with the ULIRGs in our sample, we also did 
three-blackbody fits, with both silicate and PAH emis- 
sion included, to the three average PG QSO spectra from 
Paper II. These spectra were divided by FIR strength 
into FIR-strong QSOs, FIR-weak, and FIR-undetected. 
These new fits differ from the fits in Papers I and III in 
that Paper I fit only the PAHs using simple Lorentzian 
fits to a few individual features, while Paper HI presented 
more formal fits but did not include PAH emission. Here 
we include all components, and make sure that the PAH 
fits of the QSOs are done using the same procedures ap- 
plied here to the ULIRGs. 

6. RESULTS 

The IRS spectra of all ULIRGs in the current study 
are presented in Figure [6l with archival photometry over- 
plotted: J, H, and K, and flux densities at 12, 25, 60, 
and 100 /im. The NIR photometry is from 2MASS and 
Sanders et al. (1988a). The far-infrared photometry is 
mostly from the IRAS Faint Source Catalog, but also in- 
cludes some data from Sanders et al. (2003) and a few 
ISO 12 /im points from Klaas et al. (2001). The IRS 
spectra of the QSOs were presented in Paper II, so they 
are not shown here again. The basic results from our 
analysis are fisted in Tables [51 -Tfl 

In this section, we first describe the results from our 
analysis of the broadband continuum emission in Section 
6.1 before discussing the absorption and emission- line 
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features in Section 6.2 and Sections 6.3 — 6.6, respec- 
tively. 

6.1. Broadband Continuum Emission 
6.1.1. Average Spectra 

We divided the sample into various categories, and pro- 
duced average spectra by normalizing individual spectra 
to the same rest- frame, 15 /^m flux density. 

First, we show average spectra of "cool" (/25//6O < 
0.1) ULIRGs, "warm" (/25//6O > 0.1) ULIRGs, and 
all PG QSOs (Figure [7]). The average spectrum of 
cool ULIRGs shows a steep 5-30 /zm SED with strong 
PAH features, H2 lines, and low-ionization fine struc- 
ture lines typical of starburst-dominated systems, while 
PG QSOs have a shallow 5-30 SED with silicate fea- 
tures in emission and relatively strong high-ionization 
lines and weak PAH features typical of AGN-dominated 
systems. The properties of the average spectrum of warm 
ULIRGs are intermediate between those of cool ULIRGs 
and PG QSOs. 

A similar exercise is carried out using the infrared lumi- 
nosity: \og[L{lR) / Lq] < 12.3 and > 12.3 (this threshold 
was selected to get roughly equal number of objects in 
each luminosity bin). Figure [H] shows that the contrast 
between low- and high-luminosity ULIRGs is nowhere 
near as large as between cool and warm ULIRGs. High- 
luminosity ULIRGs tend to be slightly warmer than low- 
luminosity ULIRGs so the differences we see in Figure 8 
are readily explained by the correlation between Spitzer 
spectral characteristics and /25 / /eo discussed in the pre- 
vious paragraph. 

In Figure [9l we divide the 1- Jy ULIRGs according to 
their optical spectral type. The overall 5-30 fim SED 
clearly steepens and the silicate absorption feature and 
H2 and low-ionization fine structure emission lines clearly 
become stronger as one goes from the QSOs, to the 
Seyfert Is, the Seyfert 2s, and finally to the LINER and 
H Il-like ULIRGs. The averaged spectra of these last two 
classes of ULIRGs are hardly distinguishable from each 
other with the possible exception of the silicate absorp- 
tion through, where we are S/N-limited (this is consistent 
with the /S'O-based results of Lutz et al. 1999). 

Finally, in Figures [TU] and [IH we divide the ULIRG 
sample based on the equivalent width of the PAH 7.7 
/im feature and the effective optical depth of the silicate 
absorption trough, and compare the results once again 
to the average spectrum of PG quasars. The strength 
of the PAH feature and effective optical depth of the sil- 
icate feature are derived from the SED decomposition 
described in Section 6.2 and Section 6.3. Strong absorp- 
tion features of water ice -I- hydrocarbons (5.7-7.8 /im), 
silicate (8.5 - 12 /im), C2H2 13.7 /im, and HCN 14 fun 
are detected in the absorption-dominated ULIRGs, sim- 
ilar in depth to the features seen in the heavily absorbed 
spectra of NGC 4418 and other galaxies including some 
ULIRGs (e.g., Spoon et al. 2001, 2002, 2004, 2006). Sil- 
icate absorption is visible in both PAH-dominated and 
PAH- weak systems. Similarly, PAH emission is detected 
regardless of the depth of the silicate absorption fea- 
ture. As we discuss quantitatively in Section 6.3, Section 
7.1, and Section 7.3, this lack of a clear trend between 
PAH strength and silicate absorption is largely due to 
the strong-AGN ULIRGs, which have weak PAHs and 



weak silicate absorption (see also Desai et al. 2007 and 
Spoon et al. 2007). 

6.1.2. Continuum Diagnostics 

In Figure 1121 we compare the continuum flux ratios 
/15//6, /30//6, /30//15, and /25//60 of aU ULIRGs and 
quasars in the sample. The "reddening" of the SED as 
one goes from the QSOs to the ULIRGs is evident in 
all panels of this figure. The best segregation by opti- 
cal spectral type is seen when using /25//60 and /ao/Zis- 
QSOs, Seyfert 1 ULIRGs, Seyfert 2 ULIRGs, and H H- 
like -I- LINER ULIRGs form a sequence of increasing 60- 
to-25 and 30-to-15 /im flux ratios, the H Il-like ULIRGs 
being indistinguishable from the LINER ULIRGs. Inter- 
estingly, the optically-selected starbursts observed with 
ISO (Verma et al. 2003) have 25-to-60 /im and 30-to-15 
/tm flux ratios (Brandl et al. 2006) that are intermediate 
between those of H H-like/LINER ULIRGs and Seyfert 
2 ULIRGs. We return to the /30//15 ratio in Section 
6.1.3 and in Section 7.1, where we discuss MIR spectral 
classification. 

Figure [T51 shows that PG QSOs and optically classified 
Seyfert 1, Seyfert 2, and LINER + H H-like ULIRGs 
progressively have weaker MIR emission relative to their 
FIR emission (see also Figure 7 in Paper I). The solid 
fine represents X(MIR) = L(FIR). QSOs are weU fit, on 
average, by the dotted line above it: L(MIR) w 2L(FIR). 
This line traces AGN-dominated systems and may be 
used in principle to estimate the AGN contribution to the 
ULIRG power. We return to this point in Section 6.1.3 
and Section 7.1 of this paper. Here we simply note that 
the extrapolation of this line to higher MIR luminosities 
is a good fit to the measurements of some, but not all, 
Seyfert 1 ULIRGs. These latter objects are more MIR- 
luminous than QSOs but they have only slightly cooler 
SEDs than QSOs (e.g.. Figures H and [HD . 

6.1.3. Results from SED Decomposition 

The results from the SED decomposition analysis de- 
scribed in Section 5.2 are presented in Figures [HI - I 161 
Figure [HI shows the distributions of temperatures for the 
cold, warm, and hot blackbody components used in the 
fits. Note that, as mentioned in Section 5.2, the tem- 
peratures of hot components with T > 1000 K are not 
well constrained in the fits. However, it is clear that 
Seyfert ULIRGs, particularly Seyfert 1 ULIRGs, show a 
tendency to have a warmer hot component than H Il-like 
and LINER ULIRGs. This separation is not seen in the 
warm and cold components. 

Figure [12] presents the distributions of observed 
monochromatic 8-, 15-, and 30- /im blackbody to total 
infrared luminosity ratios for all ULIRGs in the sample 
according to their optical spectral types. These black- 
body luminosities represent the sum of all blackbody 
components fitted to the IRS spectra of these objects, un- 
corrected for extinction. K-S and Kuiper tests on these 
figures confirm the stronger MIR (8 and 15 /im but not 
30 /im) continuum emission in Seyfert ULIRGs, particu- 
larly Seyfert Is, than in H Il-like or LINER ULIRGs. A 
similar result is found in Figure [TBl i. where the 5.4 - 25 
/im "pure" (PAH-free and silicate-free) blackbody emis- 
sion is compared to the FIR emission. This ratio is very 
strongly correlated with /30//15 (Figure [TBH). Table |S] 
lists for each ULIRG the observed PAH- and silicate-free 
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5-25 luminosities as a (logarithmic) fraction of the 
FIR luminosity. 

One can safely assume that the continuum emission 
from the atmospheres of the young stars in ULIRGs, and 
the very hot (10'^ K) small grain NIR dust emission com- 
ponent inferred in ISOPHOT spectra of normal galax- 
ies (Lu et al. 2003) and presumed to also exist in these 
objects, do not contribute significantly to the observed 
continuum above 5 /im. Consequently, the results in 
Figures [Ml- 1 161 most likely reflect an elevated AGN con- 
tribution to the MIR emission of Seyfert 1 and 2 ULIRGs 
relative to that of H Il-hke or LINER ULIRGs. This is 
discussed more quantitatively in Section 7.1. 

6.2. Absorption Features 

Table [HI lists the effective 9.7 fim silicate optical depth, 
r| Y, defined as 

r«^exp(-rf7) = heM-4.7l (1) 

where IcS ^ ^ h and the sum is over the blackbody 
components i. Note that the silicate feature is in emis- 
sion in four Seyfert 1 ULIRGs (Mrk 1014, F07598+6508, 
3C 273, and 21219-1757) and aU QSOs (Paper III). Also 
note that our fits were kept simple and neglected possible 
variations in the 18/10 /im absorption ratios in ULIRGs, 
so we cannot constrain the geometry of the dust dis- 
tribution in these objects in detail (e.g., Sirocky et al. 
2008; Li et al. 2008). For the same reason, we do not 
attempt to constrain the fraction of silicate absorption 
that is from crystalline silicates rather than amorphous 
silicates (e.g., Spoon et al. 2006). Finally, it is impor- 
tant to point out that the effective silicate optical depth 
discussed here is a true optical depth, defined with re- 
spect to the unextincted blackbody flux level derived 
from our fits. It is therefore different from those pub- 
lished in earlier studies (e.g., Brandl et al. 2006; Spoon 
et al. 2007; Armus et al. 2007; Imanishi et al.' 2007), 
where the depth of this feature is measured empirically 
with respect to the observed (extincted) continuum. It 
would be the same if the continuum and silicates were 
equally extincted but unfortunately that is not generally 
the case. A comparison between our measurements and 
those published in Armus et al. (10 objects) indicate that 
Tg^y ~ (2.5 ± 0.8)rg.7(Armus) with a median ratio of 2.8. 

Figure [T7| shows the distribution of Tg^ vs. the opti- 
cal spectral types of ULIRGs. The broad distribution of 
silicate strength in ULIRGs is well known from previous 
studies (e.g., Hao et al. 2007; Spoon et al. 2007). K-S and 
Kuiper tests indicate that LINER and H Il-like ULIRGs 
have significantly larger Tg^ than Seyfert ULIRGs on av- 
erage, in general agreement with Spoon et al. (2007) and 
Sirocky et al. (2008). H ILhke ULIRGs are statistically 
indistinguishable from LINER ULIRGs and the same is 
true between Seyfert 1 and Seyfert 2 ULIRGs. Combin- 
ing these results with those in Section 6.1, we find that 
all of the Seyfert 1 ULIRGs and most of the Seyfert 2 
ULIRGs cluster in the lower-left portion of the vs. 
L{8 /im)/i(IR), i(15 Mm)/£(IR), /ao/Ze, and /30//15 
diagrams (Figure [TB|) . No clear trend is seen between 
rf^ and /30//15 among HILlike and LINER ULIRGs, 
contrary to the optically-selected starburst galaxies of 
Brandl et al. (2006), where objects with strong silicate 
absorption tend to have a steeper MIR continuum. 



Also listed in Table [5] are the equivalent widths of the 
sum of the H2O ice (5.7-7.8 /im) and aliphatic hydrocar- 
bons (6.85 -|- 7.25 /im) features, and individual equivalent 
widths for C2H2 13.7 ^lm and HCN 14 ^m. [The HCO+ 
12.1 /im and HNC 21.7 /im features, whose millimetric 
transitions are important diagnostics of radiative pump- 
ing and possibly the presence of AGN (e.g., Imanishi et 
al. 2006b; Guelin et al. 2007), were not detected in any 
individual object. Upper limits of 5 x 10^^ /im and 
2 X 10"'^ /im were measured for the equivalent widths 
of IICO+ 12.1 /im and HNC 21.7 /im, respectively, in 
the average spectrum of ULIRGs with t|* > 3.86, the 
median Tg^]. The equivalent widths of C2H2 and HCN 
were measured directly, using SMART, only in objects 
with obvious detections (26 and 20 objects, respectively, 
or 35% and 27% of all ULIRGs). In contrast, the equiv- 
alent width of the H2O ice -|- hydrocarbons feature was 
derived from the SED fit of each object. This equiva- 
lent width is calculated with respect to the blackbody 
continuum only {+ silicate continuum in four Seyfert 1 
ULIRGs) i.e., PAH emission is not counted as contin- 
uum. Upper limits were set as follows. Each spectrum 
was inspected visually to determine whether or not the 
H2O absorption fit was robust. For those judged ques- 
tionable, the measured equivalent width was set as an 
upper limit. For those objects with no H2O absorption 
(often because we fixed it that way) and that have a 
significant PAH contribution (which turn out to be HII 
galaxies, LINERs, or Seyfert 2s), EW(H20 + HC) was 
assigned a limit of < 0.1 /xm. This is obviously uncer- 
tain, but implies a somewhat reasonable 5% sensitivity 
to absorption if the PAH contributes half the emission at 
these wavelengths. For the Seyfert Is with limits only, 
the upper limit was set equal to that of the lowest Seyfert 
1 measurement. Firm measurements exist for 46 objects 
(62% of all ULIRGs) and upper limits on all the others. 

Interestingly, objects with the strongest C2H2 13.7 /im 
and HCN 14 /im absorption features are not necessarily 
those with the strongest silicate and H2O ice features. 
Indeed, Figure [11] shows that the equivalent widths of 
H2O ice + hydrocarbons, C2H2, and HCN correlate only 
loosely with Tg^ and between each other. This implies 
significant variations in composition of the dense absorb- 
ing material from one ULIRG to the next. The strongest 
correlation is found between C2H2 and HCN, which a 
posteriori is not surprising since both features are be- 
lieved to be tracers of high-density (> 10* cm~'^), high- 
temperature chemistry (e.g., in Young Stellar Objects; 
Lahuis & van Dishoeck 2000; Lahuis et al. 2006, 2007). 

6.3. PAHs 

The PAH 6.2 and 7.7 /im equivalent widths and the 
total PAH to infrared and FIR luminosity ratios are listed 
in Table[71 The PAH luminosities are taken from our fits, 
and are corrected for extinction in the four sources with 
fitted PAH extinction (see Section 5.2 for more details). 
The equivalent widths are computed by dividing the PAH 
luminosity by the observed (extincted) continuum fluxes 
at 6.22 and 7.9/im. Without proper fits, the PAH 7.7/im 
equivalent width measurements are subject to errors in 
the silicate 9.7/im absorption correction. However, our 
fits to the entire 5-30/im IRS spectra take this effect into 
account in a robust manner. In what follows, we use 
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the 7.7^m feature exclusively, though the 6.2/im feature 
gives identical results. 

The results of our fits (Section 5.2) suggest that the de- 
tected PAHs in our sources are lightly extincted {Ay < 
10). This means that the extinction towards the observ- 
able PAH-emitting regions in ULIRGs is small compared 
to the sometimes heavily extincted blackbody-emitting 
regions. This is consistent with the detection of spatially 
extended PAH emission in compact U/LIRGs by Soifer 
et al. (2002). 

That said, we cannot rule out heavily extincted PAHs 
in the cores of ULIRGs. We can set limits on the contri- 
bution of such heavily extincted components to the total, 
unextincted PAH emission. First we assume that any 
heavily obscured PAH emission is extincted to the same 
degree as the continuum. Then, for the median Tgt/ 
in our sample (3.9), any heavily obscured PAH emission 
must constitute less than about a third of the total unex- 
tincted PAH emission for it to not significantly alter the 
fit. This obscured component could rise to half of the 
total unextincted emission if Tg'^/ was about twice the 
median (6 — 7) . The four sources where PAH extinction 
is detected (Section 5.2) may be cases where obscured 
PAH emission starts to dominate. 

Both PAH (e.g., Forster-Schreiber et al. 2004; Peeters 
et al. 2004; Calzetti et al. 2007) and FIR emission (e.g., 
Kennicutt 1998) are tracers of star formation in quies- 
cent and actively star forming galaxies. We also argue 
in Papers I and II that PAH and FIR emission in PG 
QSOs are produced by star formation. Wc find a fairly 
tight distribution of L(PAH)/L(FIR) in ULIRGs (Figure 
[^Ob and 6), consistent with previous studies (Peeters et 
al. 2004) as well as the notion that both trace star for- 
mation. [i(PAH) is the total PAH flux in the 5 - 30^m 
range.] K-S and Kuiper tests indicate no significant trend 
with optical spectral type. The mean and standard devi- 
ation are log i(PAH)/L(FIR) ~ -1.71 ± 0.3. The same 
conclusions apply if L{IR) is substituted for L(FIR), and 
we measure log i(PAH)/L(IR) ~ -1.96 ± 0.3. If the 
7.7/im PAH luminosity is substituted for the total lumi- 
nosity, these PAH ratios are lower by 0.4 dex. Thus, our 
results are consistent with L(PAH, 7.7/im)/L(FIR) for 
PG QSOs (-2.0 ± 0.3; Paper I). 

However, we do find that galaxies with stronger than 
average silicate absorption have smaller i(PAH) /L(FIR) 
ratios by a factor of 2 than galaxies with weaker than av- 
erage absorption, as verified with K-S and Kuiper tests 
(Figure [20h). In fact, the PAH-to-FIR ratio anticorre- 
lates with effective silicate optical depth, such that larger 
extinction corresponds to smaller L(PAH)/i(FIR) (Fig- 
ure [20li). This effect is most pronounced in the HII and 
LINER ULIRGs. As we note above, ~ half of the in- 
trinsic PAH emission may be completely buried in the 
most heavily obscured sources. A factor-of-two correc- 
tion could close at least some of the discrepancy between 
PAH-to-FIR ratios in heavily obscured and lightly ob- 
scured ULIRGs, and further tighten the distribution of 
L (PAH) /i (FIR). However, we argue below that the dif- 
ferences we observe are more likely due to a real suppres- 
sion of PAH emission. 

In Figure [2Th . we show the distribution of 7.7/im equiv- 
alent widths, which is quite broad and shows significant 
optical spectral type dependence. Seyfert 1 ULIRGs have 



much smaller PAH equivalent widths on average than 
H II ULIRGs (-0.78 vs 0.13), while the PAH equivalent 
widths of Seyfert 2 (-0.29) and LINER (-0.05) ULIRGs 
fall in between these values, confirming earlier ISO re- 
sults (e.g., Genzel et al. 1998; Lutz et al. 1999) as well 
as recent Spitzer results (e.g., Desai et al. 2007; Spoon 
et al. 2007). PG QSOs overlap with Seyferts. 

A weak luminosity dependence is also present. 
ULIRGs with log[i(IR)/iQ] > 12.4 have shghtly smaller 
7.7/im PAH equivalent widths than lower luminosity ob- 
jects (-0.36 ± 0.09 vs. 0.04 ± 0.08; Figure ETt), in 
agreement with earlier ISO results (e.g., Lutz et al. 
1998b; Tran et al. 2001). The PAH equivalent widths 
of PG QSOs are similar to those of Seyfert 1 ULIRGs, 
but they do not follow the trend with infrared lu- 
minosity of the ULIRGs (this is not surprising since 
PG QSOs were not selected through infrared methods 
like the ULIRGs). The slight IR luminosity dependence 
of EW(PAH) among ULIRGs coincides with the excess of 
Seyfert 1 ULIRGs (7 out of 9) and deficit of HII ULIRGs 
(2 out of 18) in the high-luminosity bin of our sample. In 
other words, it parallels the well-known infrared luminos- 
ity dependence of the optical spectral types of ULIRGs 
(Veilleux et al. 1995, 1999a). 

Given the trend with optical spectral type, it is not 
surprising to find that ULIRGs with warmer quasar-like 
MIR continua exhibit smaller PAH equivalent widths 
than cooler systems (Figure [2Tb). The /30//15 ratio 
is particularly efficient at separating objects, including 
PG QSOs, according to their PAH equivalent widths. 
Arguably it is even better at it than the optical spec- 
tra type, since there is a correlation between /30//15 and 
EW(PAH) among galaxies of a given spectral type. Our 
results also indicate that /30//15 is a better proxy for 
EW(PAH) than fso/fe and even fis/fe, the continuum 
color diagnostic used by Laurent et al. (2000). We return 
to this point in Section 7.1. 

As with the PAH-to-FIR ratios (Figure [Uli), the 
7.7/im equivalent width correlates strongly with extinc- 
tion in HII/LINER ULIRGs (Figure EHf) . What is the 
origin of these dependences? For the PAH-to-FIR ra- 
tio, we cannot rule out extinction effects. We can for 
EW(PAH), as long as the continuum and any unob- 
served, heavily obscured PAHs are extincted to roughly 
the same degree. The equivalent width is, however, af- 
fected by a strongly varying amount of warm continuum. 
We observe a broad distribution of 8/im-to-IR ratios (Fig- 
ure [T5li), suggesting that the 8/im continuum plays an 
important role in regulating EW(PAH). The anticorrela- 
tion of both PAH-to-FIR ratio and EW(PAH) with ex- 
tinction also points to the presence of PAH suppression 
at high extinction / low EW(PAH). This supression may 
be due to effects of high density in the cores of ULIRGs, 
or to destruction of PAHs in the harsh radiation field of 
AGN (whose importance increases with increasing opti- 
cal depth in HII/LINER ULIRGs; §7.3). In our recipe 
for computing AGN contribution from EW(PAH), we as- 
sume that (a) PAH emission is due to star formtion and 
that (b) an AGN causes both an increase in the 8/im 
continuum and PAH suppression (§7.1). 

A qualitatively similar result was found by Desai et 
al. (2007) and Spoon et al. (2007) using large samples 
of starbursts, AGN, and ULIRGs (some of these are 
also part of our sample). The HII and LINER ULIRGs 
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in our sample populate a diagonal sequence joining the 
highly absorbed, weak-PAH ULIRGs with the unab- 
sorbed, PAH-dominated systems; this sequence coincides 
with the "diagonal branch" of Spoon et al. On the other 
hand, all Seyfert 1 ULIRGs and many, but not all, Seyfert 
2 ULIRGs in our sample have both weak silicate ab- 
sorption and weak PAHs; they populate what Spoon et 
al. call the "horizontal branch". The existence of these 
two branches may reflect true intrinsic differences in the 
power source and/or nuclear dust distribution between 
galaxies on the two branches (e.g., Levenson et al. 2007; 
Spoon et al. 2007; Sirocky et al. 2008). We return to this 
point in Section 7.1 and Section 7.3. 

6.4. Fine Structure Lines 

In this section we use the strengths of the fine struc- 
ture lines to constrain the properties of the warm ionized 
gas near the central energy source of our sample galaxies. 
We first discuss the low-excitation features that are com- 
monly detected in star-forming galaxies before discussing 
the high-excitation lines, direct probes of the AGN phe- 
nomenon. 

Figure compares the luminosity of the [Ne II] 12.8 
fim line in the ULIRGs and QSOs of our sample with 
the FIR luminosity. The average and median values of 
the [Ne II] /FIR luminosity ratios are remarkably simi- 
lar regardless of the optical spectral type, including the 
Seyfert 1 ULIRGs and the QSOs: log L([Ne II])/L(FIR) 
= -3.35 ± 0.10. The similarity of this ratio for ULIRGs 
and QSOs was first pointed out in Paper I, where this 
result in combination with the similar PAH-to-FIR lumi- 
nosity ratio (§6.3) was used to argue that the bulk of the 
FIR luminosity in QSOs is produced via obscured star 
formation rather than the AGN. Not surprisingly, the 
less obscured optically-selected ISO starbursts (Verma et 
al. 2003) and Seyfert galaxies (Sturm et al. 2002) plotted 
in Figure have noticeably larger [Ne II] /FIR ratios 
(by a factor of ^ 2 and 4, respectively). 

The ([Ne III] 15.5 ^im)/{[Ne II] 12.8 /im) line ratio 
is commonly used to diagnose the excitation proper- 
ties of star-forming galaxies (e.g., Thornley et al. 2000; 
Verma et al. 2003; Brandl et al. 2006). Since the ion- 
ization potentials needed to produce Ne"*" and Ne"*""*" are 
21.6 and 41.0 eV, respectively, the [Ne III]/[Ne II] ra- 
tio is sensitive to the hardness of the ionizing radia- 
tion and therefore to the (effective temperature of the) 
most massive stars in a starburst or to the presence of 
an AGN, if applicable. Figure [53] shows this ratio as a 
function of the infrared and FIR luminosities and the 
MIR continuum colors, /25//6O and /30//15. No ob- 
vious trend with the F/IR luminosities is seen among 
ULIRGs. However, a clear dependence is seen with opti- 
cal spectral type and MIR continuum colors, confirm- 
ing earlier studies (e.g.. Dale et al. 2006; Farrah et 
al. 2007). Larger [Ne III]/[Ne II] ratios go hand- in- 
hand with warmer MIR continuum. QSOs have larger 
[Ne III]/[Ne II] ratios on average than Seyfert ULIRGs, 
and Seyfert ULIRGs have larger ratios on average than 
HE-like and LINER ULIRGs. K-S and Kuiper tests in- 
dicate that the [Ne III]/[Ne II] ratios of LINER ULIRGs 
are statistically indistinguishable from those of Hll-like 
ULIRGs, and the same statement also applies when com- 
paring Seyfert 1 ULIRGs with Seyfert 2 ULIRGs. The 
two Seyfert 1 ULIRGs with [Ne III] upper limits are 



F07598-I-6508 and F13218-H0552. Both of them have un- 
usually small optical narrow-line [OIII] A5007/H/3 ratios, 
consistent with low-luminosity high-excitation regions 
(e.g., Kim et al. 1998). The dependence of [Ne III]/[Ne II] 
on spectral type and MIR colors induces a slight trend 
between this ratio and EW(PAH 7.7) (Figure [23;). No 
obvious trend between [Ne III]/[Ne II] and EW(PAH 7.7) 
is seen within Hll-like and LINER ULIRGs, a result that 
is consistent with that found for optically-selected star- 
burst galaxies (Brandl et al. 2006). Given the relatively 
high metallicity of ULIRGs (Section 6.6 in this paper and 
Rupke et al. 2008), the large [Ne III]/[Ne II] ratios seen 
in Seyfert ULIRGs cannot be explained by star formation 
alone as in the case of low-metallicity dwarf galaxies. 

Other low ionization fine structure lines such as [Fe II] 
25.99 fj.m, [S III] 33.48 ^im, and [Si II] 34.82 fj.m have 
been found to be useful diagnostics of activity in galactic 
nuclei (e.g., Lutz et al. 2003; Sturm et al. 2005; Dale et 
al. 2006). Unfortunately, these lines are often redshifted 
out of the wavelength range of our data so they cannot 
be used for any kind of statistical analysis. We do not 
discuss these lines any further in this paper. 

The ionizing spectra of all but the hottest O stars cut 
off near the He II edge (54.4 eV), so the detection of 
[O IV] 25.9 /im from three-times ionized oxygen with 
ionization energy % ~ 55 eV is potentially a good in- 
dicator of AGN activity. This line is detected in 30/34 
QSOs, 3/9 Seyfert 1 ULIRGs, 8/13 Seyfert 2 ULIRGs, 
6/28 LINER ULIRGs, and only one of the 18 HH-like 
ULIRGs (F21208-0519:N) with high-resolution spectra. 
In Figure we plot the (upper hmits on the) [O IV] 
25.9 /im luminosity versus the FIR luminosity of ULIRGs 
and QSOs. The solid line is a fit to the data of HII 
ULIRGs, so it is formally only an upper limit, as indi- 
cated by the arrows. This upper limit is above, therefore 
consistent with, the measured values in ISO starbursts. 
The PG QSOs and the few Seyfert 1 ULIRGs with [O IV] 
detections lie on average ~1.2 dex above that line. All 
Seyfert 2 ULIRGs with [O IV] detections he ~ 0.8 dex 
above that line, while the upper limits on [O IV] de- 
rived for the other Seyfert ULIRGs are consistent with 
those for the HII ULIRGs and reflect the flux detection 
threshold across the sample. Interestingly, the optically- 
selected ISO Seyfert galaxies have [O IV] /FIR luminosity 
ratios that are similar to those of the PG QSOs. 

As in the case of [Ne III]/[Ne II], there is no strong 
trend between the [O IV] /[Ne II] ratios and M/IR lumi- 
nosities of ULIRGs, but a strong dependence with optical 
spectral type and MIR continuum colors is detected (Fig- 
ure [M]). These results are similar to those found by Far- 
rah et al. (2007) on a different but overlapping sample of 
ULIRGs. The lack of an obvious luminosity dependence 
among ULIRGs may be surprising in the light of the 
optical and ISO results which suggest a larger AGN con- 
tribution to the bolometric luminosity among the more 
luminous ULIRGs (e.g., Veilleux et al. 1995, 1999a; Lutz 
et al. 1998b; Tran et al. 2001). Possible explanations for 
this apparent discrepancy include: (1) the optical spec- 
tral classification is affected by dust obscuration so it is 
not reliable; (2) [O IV] 25.9 fim is not as good an AGN 
indicator as the PAH equivalent width (e.g., contaminat- 
ing [O IV] emission from WR stars and ionizing shocks, 
Lutz et al. 1998a; Abel & Satyapal 2008); (3) given typ- 
ical ULIRG redshifts and actual IRS sensitivity at > 30 
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/iin, detecting [O IV] is difficult. The number of ULIRGs 
with actual [O IV] detection is small, especially among 
Hll-like and LINER ULIRGs, so small number statistics 
mask the correlation. We favor this last possibility. Ex- 
planation #1 can be rejected outright since we detect in 
the present paper (as in Lutz et al. 1999) clear correla- 
tions of the MIR parameters with optical spectral types. 
If the optical classification were unreliable, these rela- 
tions would be erased. Scenario #2 seems unlikely since 
starbursts with large [O IV] /[Ne II] are (low metallicity) 
dwarfs; at the relatively high metallicity of ULIRGs, this 
ratio is expected to be small. Note, however, that sev- 
eral Seyfert 1 ULIRGs have no detected [O IV] emission. 
This is not purely a sensitivity effect. As pointed out in 
our discussion of the [Ne III]/[Ne II] ratios in Seyfert 1 
ULIRGs, there is strong corraborating evidence at op- 
tical wavelengths that many of these objects have high- 
excitation regions of unexpectedly low luminosity. The 
exact cause of this effect is unclear. With these caveats 
in mind, we will use the [O IV] /[Ne II] ratio in Section 
7.1 as a diagnostic of nuclear activity in ULIRGs. 

Despite the fact that [Ne V] 14.3/im is fainter in these 
sources than [O IV] , the better S /N ratio of the IRS data 
at shorter wavelengths has allowed us to put similar con- 
straints on the [Ne V] and [O IV] lines. ([Ne V] 24.3/L(m 
was also detected in many sources with [Ne V] 14.3/L(m 
emission, but at a lower rate than [Ne V] 14.3 /im). 
[Ne V] lA.Sfim was detected in 25/34 QSOs, 4/9 Seyfert 1 
ULIRGs, 9/13 Seyfert 2 ULIRGs, 4/28 LINER ULIRGs 
(F04103-2838, UGC 5101, F13335-2612, NGC 6240), 
and even one of the 18 HII ULIRGs with high-resolution 
spectra (F20414-1651). The relatively small redshifts of 
three of the five detected HII/LINER ULIRGs makes it 
apparent that sensitivity plays a role in the detectabil- 
ity of these lines in sources where they are intrinsically 
weak. Once again, the relatively modest number of de- 
tections in Seyfert 1 ULIRGs point to intrinsically weak 
high-excitation regions in some of these objects. 

The very high ionization potential of Ne*+, x = 97.1 
eV, makes this line an unambiguous signature of nuclear 
activity. Indeed the separation with spectral type and 
MIR continuum colors previously seen in [O IV]/(F/IR) 
and [O IV] /[Ne II] is clearer when [Ne V] is sulDstituted 
for [O IV] (Figure and [25]). This separation is further 
emphasized in Figure [26l where we compare the values of 
[Ne V]/[Ne II] with [O IV] /[No II] measured in our sample 
of QSOs and ULIRGs. The solid diagonal line in these 
diagrams is a line of constant [Ne V]/[0 IV] and changing 
Ne II]. This line may be interpreted as a mixing line if 
Ne V and [O IV] are only produced by an AGN and 
Ne II] by starburst activity. The tickmarks along the 
fine indicate the percentage contribution of the starburst 
to [Ne II] (from to 99%). Optically-selected ISO Seyfert 
galaxies lie in the same region as the PG QSOs in all 
these diagrams (Figures UHl and . We will return 
to this last figure in our discussion of the energy source 
in ULIRGs and QSOs (Section 7.1). 

In the ULIRGs and PG QSOs with high-S/N SL spec- 
tra, we also searched for redshiftcd [No VI] 7.65 /im (x = 
126 eV), another powerful AGN indicator. This line was 
unambiguously detected in 6 QSOs, 2 Seyfert 1 ULIRGs, 
and 5 Seyfert 2 ULIRGs, but in none of the LINER and 
Hll-like ULIRGs. Confusion between [Ne VI] and PAH 
substructure can mask weak [Ne VI] emission in the lat- 



ter objects. By and large, [Ne VI] follows the same trends 
with spectral type and MIR continuum colors as [Ne V] 
(Figures [2li and[27l). 

The [Ne VI] /[O IV] vs. [Ne VI] /[Ne II] diagram from 
Sturm et al. (2002) is reproduced in Figure [^5b . Both 
axes scale with the AGN excitation, but [Ne VI] /[Ne II] 
can also be influenced by contributions from star forming 
regions to the [Ne II] line. As shown by the grid of AGN 
models from Groves et al. (2004) , pure AGN are expected 
to lie roughly along a diagonal line in this diagram (both 
ratios increase with increasing hardness of the radiation) . 
Composite sources, however, have stronger [Ne II] lines 
than pure AGN, so they are expected to lie to the left 
of the pure AGN sources. Because Seyfert ULIRGs are 
composite sources, with significant starburst contribu- 
tion to [Ne II] (Figure [55]) , they lie leftward of most of 
the comparison Seyferts and PG QSOs. The compari- 
son Seyferts are in decent agreement with the models, 
though some may suffer minor starburst contamination 
to [Ne II]. The [Ne VI] /[Ne II] ratios of most of the 
PG QSOs are too high by factors of 2 — 3. There is some 
disagreement in [Ne VI] /[O IV], as well, which may be 
due to incorrect Ne/0 abundance ratios in the models 
(§6.6). 

A much better match with the models is found when 
considering the [Ne VI] /[Ne III] vs. [Ne VI] /[Ne V] di- 
agnostic diagram (Figure [28b ). which is independent of 
relative metal abundance effects. This suggests that the 
bulk of the [Ne VI], [Ne V], and possibly even [Ne III] 
emission in PG QSOs, Seyfert ULIRGs, and ISO Seyferts 
is produced by the AGN. 

6.5. Molecular Hydrogen Lines 

Three lines from rotational transitions of warm II2 were 
regularly detected in the spectra of ULIRGs and QSOs: 
w = 0./== 3-1 S(l) 17.04 ^m, i; = J = 4-2S(2) 12.28 
/xm, and w = J = 5 - 3 S(3) 9.67 /im. The w = J = 
2 — S(0) 28.22 /xm transition was detected in 10 objects 
(F09039-f 0503, UGC 5101, F12112+0305, F13335-2612, 
Mrk 273, F14248-1447, F21208-0519:N, PG 1211+143, 
PG 1440+356, and B2 2201+31 A). Higher-level tran- 
sitions were also detected in a few objects [S(4), S(5), 
S(6), and S(7) in F09039+0503 and F15130-1958, and 
S(5) in UGC 5101, F12112+0305, F17208-0014, and 
PG 1700+518]. 

To first order, the strengths of the H2 lines scale lin- 
early with the star formation rate indicators of ULIRGs: 
F/IR, [Ne II], and PAH luminosities (Figure [M]) . This is 
true in detail for the Hll-like ULIRGs, but strong depar- 
tures from the linear relation are seen among the other 
ULIRGs and the QSOs. The H2 line emission in these 
latter objects tends to be overluminous for a given F/IR 
or PAH luminosity. Similar departures from the lin- 
ear relation were seen among the SINGS LINER/Seyfert 
targets (Roussel et al. 2007) and attributed to shock 
heating. The fact that the departures are strongest 
among the QSOs of our sample suggests that heating 
by the AGN is important in these objects. Rigopoulou 
et al. (2002) came to a similar conclusion based on the 
H2 /PAH ratio of a sample of 9 optically-selected Seyfert 
galaxies. 

The temperature-sensitive H2 S(2)/S(l), S(3)/S(l), 
and S(3)/S(2) ratios may be used to shed some light 
on the possible role of the AGN. Their distributions 
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are shown in Figure [3D1 Warm Seyfert ULIRGs tend 
to have larger [smaller] S(3)/S(2) [S(2)/S(l)] ratios than 
cool HII-like/LINER ULIRGs, while no statistically sig- 
nificant trend is seen in the S(3)/S(l) ratio distribution. 
The K-S and Kuiper probabilities that the distributions 
of S(3)/S(2) and S(2)/S(l) ratios, when put in two bins 
above and below the median /25//6OJ arise from the same 
parent distribution are <2%, confirming the apparent 
trends. The number of QSOs with reliable H2 line ratios 
is too small to be able to detect statistically significant 
trends. 

We used these ratios, when available, to construct an 
excitation diagram for each object, assuming LTE, an 
ortho-to-para ratio of 3, and no extinction. Figure [31] 
shows a few examples. A straight line in these diagrams 
indicates that a single temperature applies to all tran- 
sitions, with the excitation temperature (Tex) being the 
reciprocal of the slope. For an ortho-to-para ratio of 3, 
temperatures derived from adjacent lines should increase 
with J. This is illustrated in Figure [32l where we com- 
pare Tex (J = 4-3) and Tex{J = 5-4), the excitation 
temperatures derived from the S(2)/S(l) and S(3)/S(2) 
ratios, respectively, as a function of /25//6O and the ef- 
fective silicate optical depth, A temperature 
difference that is negative implies that extinction and/or 
ortho-para effects are at play, as illustrated by the down- 
ward arrows on the right in each diagram (see caption to 
this figure for an explanation of the arrows). 

The small number of Seyfert Is and QSOs for which 
both TexiJ = 4-3) and TgxiJ = 5-4) are available pre- 
vents us from detecting any obvious dependence on op- 
tical spectral type. A visual inspection of Figure [32b 
suggests a possible trend of increasing Tex{J = 5 - 4) - 
Tex{3 = 4-3) with increasing /25//6O1 which would sup- 
port the role of AGN in heating the molecular gas in these 
objects. In fact, the K-S and Kuiper probabilities that 
the distributions of temperature differences binned ac- 
cording to /25//6O (above and below the median /25//60) 
arise from the same parent distribution are 0.6% and 6%, 
respectively, so there is a statistically significant trend. 
However, much of this trend may be due primarily to 
extinction, as shown in Figure [5^ . The K-S and Kuiper 
probabilities that the distributions of temperature differ- 
ences binned according to 7-cffcctivc (above and below the 
median Tgy'"^*'™) arise from the same parent distribution 
are only 0.8% and 7%, respectively. Thus, an impor- 
tant conclusion of this discussion is that any trends with 
/25//6O could be masked by extinction and ortho-para 
effects (c.f. Higdon et al. 2006). 

Table [8] lists the warm II2 masses derived from the 
strength of the S(l) line and the average T^x for each 
object. Values range from 0.5 to 20 x 10^ Mq with an 
average (median) of - 3.8 (3.3) and 3.6 (3.2) x 10^ Mq 
for the ULIRGs and QSOs, respectively. These values 
are slightly larger on average than those of Higdon et al. 
(2006) and imply that the warm gas mass is typically 
a few percent of the cold gas mass derived from ^^CO 
observations (0.4 - 1.5 x 10^° Mq; e.g., Solomon et al. 
1997; Downes & Solomon 1998; Evans et al. 2001, 2002, 
2006; ScoviUe et al. 2003). 



6.6. Metal Abundance 



In principle, one can use the strengths of the fine struc- 
ture lines relative to the hydrogen recombination lines in 
our data to derive the metallicity of the gas producing 
these emission features. In practice, the only hydrogen 
line within the wavelength range of our data is the very 
faint Hua 12.4 /im (H 7-6), so the S/N of our data only 
allow us to marginally detect or put upper limits on the 
strength of this line in individual objects. However, Hua 
is detected at S/N~ 6 in the average spectrum of 27 
PAH-dominated ULIRGs (Figure 36a), so we can use 
the strength of this line to derive an average metallicity 
in these systems. We follow the methods of Verma et al. 
(2003), using the ratios [Ne II] 12.8 /xm/Hua = 62 and 
[Ne HI] 15.5 /Ltm/Hua = 18 measured from the average 
spectrum to derive the abundance of neon in these sys- 
tems (Figure 366). The recombination line Hua, [Ne II] 
12.8 fim which is tracing the dominant singly ionized 
state of neon, and [Ne III] 15.5 /im tracing doubly ion- 
ized neon are found at similar MIR wavelengths. They 
can be used to reach optically obscured regions and ob- 
tain a metallicity measurement that is much less sensi- 
tive to extinction effects than results obtained in combi- 
nation of MIR lines with NIR recombination lines. We 
did not apply an extinction correction to the observed 
MIR line ratios. Adopting an electron temperature of 
5000 K appropriate for dusty starbursts (e.g. Puxley et 
al. 1989), Hua emissivity from Storey & Hummer (1995) 
and neon collision strengths from Saraph & TuUy (1994) 
and Butler & Zcippcn (1994), we find a neon abundance 
12+log(Ne/H) = 8.30. As discussed in the last para- 
graph of the present section, the value of the solar neon 
abundance is currently the subject of a heated debate. If 
we adopt the revised solar photospheric neon abundance 
of Asplund et al. (2004), our neon abundance is ~ 2.9 x 
solar. 

An underabundance compared to local luminosity- 
metallicity and mass-metallicity relations of galaxies was 
recently reported in the optical study of 100 star-forming 
LIRGs and ULIRGs by Rupke et al. (2008), who at- 
tributed it to a combination of two effects: a decrease 
of abundance with increasing radius in the progenitor 
galaxies and strong, interaction- or merger-induced gas 
inflow into the galaxy nucleus. Thirteen objects from 
the sample of Rupke et al. (2008) are in common with 
the current Spitzer sample. The oxygen abundance, 12 
+ log(0/H), of the ULIRGs in common with both sam- 
ples ranges from 8.43 to 9.04. Using the value from As- 
plund et al. (2004) for the solar oxygen abundance, 12 -|- 
log(0/H)Q = 8.66, these numbers translate into 0.6 - 2.4 
X solar. Given this relatively narrow range of abundance 
and small number of objects, it is perhaps not surprising 
that no trend was found in our sample between the oxy- 
gen abundances of Rupke et al. (2008) and any Spitzer- 
derived continuum or line ratios. In particular, we note 
that the oxygen abundance of these ULIRGs is well above 
the threshold abundance, 12 -I- log(0/H) = 8.1 or ~ 
0.3 solar, below which PAH emission is apparently sup- 
pressed (e.g., Engelbracht et al. 2005; O'Halloran et al. 
2006; Smith et al. 2007). 

Supersolar neon abundance derived from the MIR 
spectra (if the Asplund et al. 2004 value of the solar neon 
abundance is correct) and close to solar oxygen abun- 
dance derived from the optical spectra may trace differ- 
ent layers of the ULIRGs, in both extinction and abun- 
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dance. In the picture outlined by Rupke et al. (2008), it 
is plausible that less obscured regions are dominated by 
lower metallicity gas transported in from the outskirts of 
the galaxies. In contrast, the dusty inner regions may 
be dominated by more pre-enriched gas from the inner 
regions of the progenitor galaxies, compressed to the im- 
mediate circumnuclear region during the merger process 
and enriched further by the intense circumnuclear star 
formation. However, the excellent overall agreement re- 
ported in Section 7.1 between optical and MIR diagnos- 
tics of nuclear activity do not seem to favor this picture. 
These results imply that the optical line spectrum in 
most cases traces gas that "knows" what the true power 
source is and therefore should also trace gas that fairly 
samples the metallicity. An alternative explanation for 
the higher neon abundance - and the one we favor - is 
that it reflects in-situ enrichment in the most heavily ob- 
scured (densest) star-forming regions, but these regions 
are distributed throughout the ULIRG rather than pref- 
erentially near the center. 

An important caveat in comparing the optical and MIR 
metallicity measurements is the assumed solar Ne/0 ra- 
tio. The exact value of this ratio has been the subject 
of a heated debate in recent years, some groups arguing 
that it is considerably higher than the standard value 
(e.g., Drake & Testa 2005; Wang & Liu 2008, although 
see Schmelz et al. 2005 for a counterexample). A larger 
solar Nc/0 ratio would bring our Spitzer measurements 
in closer agreement with the optical results. 

7. DISCUSSION 

7.1. Energy Source: Starburst vs. AGN 

In this section, we use the data presented in Section 
6 to estimate the fractional contribution of nuclear ac- 
tivity to the bolometric luminosity of the ULIRGs and 
PG QSOs in our sample (herafter called the "AGN con- 
tribution" for short). We use six different methods based 
on (1) the [O IV] 25.9 ^m/[Ne II] 12.8 ^im ratio, (2) the 
[Ne V] 14.3 ^m/[Ne II] 12.8 /im ratio, (3) the equivalent 
width of PAH 7.7 ^m, (4) the PAH (5.9 - 6.8 ^m) to 
continuum (5.1 - 6.8 /xm) flux ratio combined with the 
continuum (14 - 15 pan) / (5.1 - 5.8 /im) flux ratio (see 
Figure 34), (5) the MIR blackbody to FIR flux ratio, and 
(6) the /30//15 continuum flux ratio. These methods are 
described in detail in Appendix A. The zero points, bolo- 
metric corrections, and basic results from each method 
are listed in Tables [9]- l 121 We compare the results from 
the various methods and look for trends with optical and 
infrared parameters in Section 7.1.1 and Section 7.1.2, 
respectively. 

7.1.1. Comparisons of Results from Different Methods 

Tables [11] H 121 and Figure [35] indicate a remarkably 
good agreement between the AGN fractional contribu- 
tions to the bolometric luminosities of ULIRGs and 
PG QSOs derived from the various methods. The mean 
ULIRG AGN contribution is ^ 38.8 ± 21.1% averaged 
over all ULIRGs and all methods (in Table 1111 the 
average-of-averages and standard errors are calculated 
by first averaging over all methods for individual objects, 
then averaging over objects). This mean ULIRG AGN 
contribution is in agreement with, and refines the results 
of, Genzel et al. (1998): ULIRGs are composite objects, 



hut on average powered mostly by star formation. The 
various methods give average AGN contributions that are 
within ~ ±10—15% of each other, taking into account the 
range of AGN fractional contributions derived from the 
fine structure line ratios measured from average spectra 
(see discussion in Appendix A). These small differences 
between the various methods can easily be explained by 
uncertainties on the pure-starburst zero points (see dis- 
cussion in Appendix A; the pure-AGN zero points are 
considered more robust since they are based on the FIR- 
undetected PG QSOs) and modest differential extinction 
[Ay ^ 10 mags) between the inner line-emitting region 
(where the bulk of the [O IV] and [Ne V] emission is pro- 
duced on average) and outer line-emitting region (where 
the bulk of the [Ne II] and PAH emission is produced 
on average). The good agreement between the various 
methods is not in contradiction with the results of Ar- 
mus et al. (2007) since here we compare AGN fractional 
contributions to the bolometric luminosities, while Ar- 
mus et al. did not apply bolometric corrections to their 
numbers so they were comparing AGN fractional contri- 
butions to the [Ne II] and MIR luminosities and found 
them to be different. 

Note that there is systematic uncertainty associated 
with the choice of what defines an AGN or starburst (or 
HII region / PDR in the case of the Laurent method; 
see Appendix A for detailed discussion on the choices 
of zero points and bolometric correction factors for each 
method). For instance, there may be a range of possible 
emission-line ratios or continua that define a "pure" AGN 
or starburst. Experiments show that reasonable changes 
in zero-point values do not significantly change the re- 
sults for a given method. Nonetheless, this uncertainty 
may contribute to the scatter observed when comparing 
differing diagnostics. To smooth over these possible sys- 
tematics, in what follows we compute the average AGN 
contribution over all methods for each object. This mini- 
mizes the chance that a stronger systematic effect in any 
one method will affect the results. 

Furthermore, we cannot rule out the possibility of a 
third class of physically distinct systems. In other words, 
a pure starburst or pure AGN may not describe all of pa- 
rameter space. For instance, one possibility is that heav- 
ily obscured systems host unique physical conditions in 
high-density cores that do not replicate starbursting or 
AGN ULIRG environments. To examine this possibility, 
we highlight systems with effective silicate optical depths 
above the median in Figure 1351 It is evident from this 
figure that HII and LINER ULIRGs with higher obscu- 
ration tend to have higher AGN contribution. Whether 
this is due to fundamentally different physics, or simply 
an obscured AGN, is unclear from this diagram. We re- 
turn to this issue in Section 7.1.2 and Section 7.3, where 
we uncover smooth trends between AGN contribution, 
obscuration, and merger phase which are difficult to ex- 
plain if fundamentally different physics were at play. 

Finally, we cannot formally rule out the possibility of 
deeply buried AGN invisible at MIR wavelengths but 
contributing significantly to the FIR emission in some 
of these objects. However, it is now considered a highly 
contrived scenario given the good agreement between the 
variety of methods used to evaluate the AGN contribu- 
tion to the bolometric luminosity. As described in Ap- 
pendix A, these methods use the full gamut of diagnostic 
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tools available at 6 - 30 um. The diagnostic features are 
produced under different conditions (density, dust con- 
tent) and over a range of distances from the center. They 
also cover a broad range in wavelength and therefore dust 
optical depth. If obscured AGN are contributing signifi- 
cantly to the FIR emission of several of these sources, one 
would expect diagnostics that use long-wavelength emis- 
sion and probe deep into the cores (e.g., /30//15 ratio) 
to give systematically different results than the others. 
This is not seen in our data. 

7.1.2. Trends with Optical Spectral Types, /25//60 ratios, 
Infrared Luminosities, and Extinctions 

We detect strong correlations between Spzteer-derived 
AGN contributions on the one hand and optical spectral 
types and /25//60 ratios on the other (Figure [36h and c). 
These results confirm and expand on earlier results. The 
AGN contribution ranges from ^^15— 35% among HII and 
LINER ULIRGs (taking into account the range of AGN 
fractional contributions derived from the fine structure 
line ratios measured from average spectra [see discussion 
in Appendix A]) to ~ 50 and 75% among Seyfert 2 and 
Seyfert 1 ULIRGs, respectively. The presence of a dom- 
inant AGN in Seyfert 1 ULIRGs was first deduced from 
the strengths of the optical/NIR broad lines in a few ob- 
jects (Figures 4 and 5 of Veilleux et al. 1997 and 1999b, 
respectively) ; the new Spitzer results now show that this 
statement applies to Seyfert 1 ULIRGs in general. The 
excellent correlation between optical spectral types and 
7.7 /im PAH-derived AGN contribution was first pointed 
out by Taniguchi et al. (1999) and Lutz et al. (1999) us- 
ing ISO data, but we have now quantified this correlation 
and detected similar ones when using the fine-structure 
line and continuum slope methods. The correlation be- 
tween AGN contributions and /25//60 ratios is equally 
strong and quantitatively confirms the qualitative state- 
ment made more than twenty years ago by de Grijp et 
al. (1985) that this ratio is an excellent indicator of AGN 
activity. The AGN contribution among cool ULIRGs 
(/25//6O < 0.2) is ~ 30% on average compared with 
60% among warm ULIRGs. 

Figure [32] also displays the AGN contributions of the 
PG QSOs. These fall right along the extrapolation of 
the spectral type and /25//60 sequences, with AGN con- 
tributions typically larger than ~ 80% among the QSOs. 
(Recall that only 8 PG QSOs - only those that are FIR- 
undetected - were used to set the pure- AGN zero points 
so this last statement is not circular.) These results bring 
support to the concept of an excitation sequence be- 
tween the cool, HII/LINER ULIRGs, the warm Seyfert- 
like ULIRGs, and the PG QSOs. They are also consis- 
tent with the evolution scenario proposed by Sanders et 
al. (1988a, 1988b), if the excitation sequence is also a 
merger sequence. This question is examined in Section 
7.3 below. 

A weaker correlation is seen between the AGN con- 
tributions and infrared luminosities of ULIRGs (Figure 
[36b). We observe average AGN contributions of ~ 34% 
and - 48% for ULIRGs with log[L(IR)/L0] below and 
above 12.4, respectively. This general trend with in- 
frared luminosity is consistent with the optical results 
of Veilleux et al. (1995, 1999a) and the ISO results of 
Lutz et al. (1999) and Tran et al. (2001). The PG QSOs 
are distinctly less infrared luminous than the ULIRGs, 



yet they have larger AGN contributions. If the evolution 
scenario of Sanders et al. (1988a, 1988b) is to apply to 
PG QSOs, the infrared-luminous starburst in these ob- 
jects must have subsided from its peak activity during 
the ULIRG phase. 

Another prediction of this evolutionary scenario is that 
the AGN eventually emerges out of its dusty cocoon. 
The only diagnostic tool at our disposal to estimate the 
amount of dust in these systems is the effective silicate 
optical depth (Section 6.2). We return to Figure [?IH . 
this time considering the AGN contribution rather than 
simply the optical spectral types. The results are shown 
in Figure [57b and summarized in Table [T5I We find a re- 
markably strong trend in AGN contribution, leading from 
the lower right through the upper region and ending in the 
lower left (we have labeled these regions RI, R2, and R3 
for convenience) . All of the objects in Rl are starburst- 
dominated. In R2, the objects have larger AGN contri- 
butions, but are still mostly starburst- dominated. In R3, 
the objects are either A GN dominated or show a balance 
between starburst and A GN. These results are consistent 
with the evolution scenario if the objects on the Spoon et 
al. diagonal branch (regions Rl and R2) are in an earlier 
phase of ULIRG evolution than objects on the left tip of 
the horizontal branch (region R3). Differences between 
ULIRGs populating Rl and R2 may also be explained 
in the context of the evolution scenario if extinction in- 
creases during the intermediate stages of merger evolu- 
tion (from Rl to R2) before dust is destroyed or blown 
away by the AGN (R3). We explore this possibility in 
Section 1.3. 

7.2. Black Hole Growth Rate 

Here we calculate the Eddington ratio, rj, i.e., the ra- 
tio of AGN bolometric luminosity to the Eddington lu- 
minosity, L(Edd) = 3.3 X 10"* (Mbh/Mq) Lq, for each 
system. The results are shown in Tables [T3| - I 151 and 
Figure [211 Two methods were used to estimate the black 
hole masses in these systems: (1) "dynamical" black hole 
masses based on the stellar velocity dispersion of the 
spheroidal component in these objects from Dasyra et 
al. (2006a, 2006b, 2007) and the stellar velocity disper- 
sion — black hole mass relation of Tremaine et al. (2002), 
(2) "photometric" black hole masses based on measure- 
ments of the H-band luminosity of the spheroidal com- 
ponent in these systems (free of the central point source) 
from Veilleux et al. (2002, 2006, 2009) and the H-band 
spheroid luminosity — black hole mass relation of Mar- 
coni & Hunt (2003). 

Photometric black hole masses are available for all 
ULIRGs and PG QSOs in the Spitzer sample, while dy- 
namical estimates are available for only a third of the 
sample. Note also that the dynamical black hole mass 
measurements of ULIRGs and PG QSOs are smaller on 
average than the photometric estimates, hence the Ed- 
dington ratios derived from the dynamical black hole 
masses are larger on average than those based on the pho- 
tometric method: for the ULIRGs in our sample, log(77) 
= -1.08 ± 0.40 and -0.35 ± 0.63 (Tables [Ttj -Qll in 
these tables, the average-of-averages and standard errors 
are calculated by first averaging over all methods for in- 
dividual objects, then averaging over objects). A similar 
discrepancy is found among the PG QSOs when com- 
paring the dynamical estimates with those from rever- 
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beration mapping or the virial method. A detailed com- 
parison of the various black mass estimates in ULIRGs 
and PG QSOs is beyond the scope of the present paper 
(interested readers should refer to Veilleux et al. 2009 
for a more detailed discussion and a table of the black 
hole masses from the various methods). Suffice it to say 
that the absolute values of all Eddington ratios are quite 
uncertain so the present discussion focuses on overall rel- 
ative trends, which should be much more robust. 

Figure [39] shows the distribution of the 
photometrically-derived Eddington ratios as a function 
of spectral types, infrared luminosities, /25//6O ratios, 
nuclear separations, and interaction classes (Table 
1141 summarizes the results). No obvious correlations 
exist between the Eddington ratios and any of these 
parameters. Weak trends may be present with the 
morphological quantities: Eddington ratio appears to 
be larger at the smallest nuclear separations and latest 
interaction classes (see Section 7.3 for description of 
interaction classes). In both cases, the addition of 
PG QSOs seems to either extend or enhance these 
trends. However, a rigorous statistical analysis of these 
data cannot confirm the trends involving the ULIRGs. 
Similarly, the number of objects with dynamical Edding- 
ton ratios (Table [T5)l is generally too small to allow us 
to detect any significant trends involving this quantity. 

7.3. Merger Evolution 

Virtually all 1-Jy ULIRGs and most PG QSOs show 
clear signs of strong tidal interaction/merger. ULIRGs 
are on-going mergers that sample the Toomre merger 
sequence beyond the first peri-passage (Veilleux et al. 
2002, 2006), while many PG QSOs are advanced mergers 
where the nuclei of the progenitor galaxies have appar- 
ently coalesced (e.g., Surace et al. 2001; Guyon et al. 
2006; Veilleux et al. 2009). It is therefore natural to ask 
whether the excitation sequence we see in Figure ISHl in 
fact corresponds to the final stages of the evolution se- 
quence first suggested by Sanders et al. (1988a, 1988b). 
In this section we examine this question by using a num- 
ber of morphological indicators of merger phase: appar- 
ent (projected) nuclear separation, lengths of tidal tails, 
compactness of merger remnant, and strength of tidally- 
induced morphological anamolies in coalesced systems. 
We also use the interaction classes of 1-Jy ULIRGs and 
PG QSOs derived by Veilleux et al. (2002, 2006, 2009) 
and based on the classification scheme of Surace (1998). 
This scheme combines all morphological indicators of 
merger phase and compares the results with published 
numerical simulations of mergers (which we describe in 
more detail below). In brief. Classes I through V corre- 
spond to first approach, first contact, pre-merger [sub- 
divided into a and b for wide (> 10 kpc) and close (< 
10 kpc) pairs], merger (subdivided into a and b for dif- 
fuse and compact systems), and old merger, respectively. 
The results are summarized in Table [TT] and presented in 
Figure [36li and e. 

A simple comparison between binary- and single- 
nucleus ULIRGs reveals a slight difference between their 
respective AGN fractional contribution to the bolometric 
luminosity (32% for the binaries versus 46% for the sin- 
gles; Table [TT|). A closer look at this result indicates that 
the increase in AGN contribution generally takes place 
when the apparent nuclear separation is less than '--^l kpc 



(Figure I36| j) . This trend with nuclear separation seems 
driven primarily by the large number of Seyfert 1 ULIRG 
among merged systems. Projection effects undoubtedly 
add scatter to the data. This is also illustrated in Figure 
136b . where we substituted the more physically meaning- 
ful interaction class for the apparent nuclear separation. 
By and large, the results on the ULIRGs are consistent 
with the evolution scenario of an increasingly more dom- 
inant AGN among late mergers, although with consider- 
able scatter. Some of this scatter may be due to multiple 
episodes of AGN dominance during the merger process. 
Most wide {NS > 6 kpc) binaries do not have confirmed 
redshifts for both pair members so the slightly larger 
AGN contributions among this class of objects may be 
due to misidentifications. Our results on the PG QSOs 
bring additional support to the evolution scenario: the 
high AGN contribution and unresolved morphologies of 
the PG QSOS falls right along the trends observed among 
ULIRGs (Figure EIH and e). 

Figure [59H and e suggest that pre-merger ULIRGs are 
indeed on average less actively accreting matter onto the 
black holes than late mergers, based on their Edding- 
ton ratios, with the PG QSOs nicely falling along these 
trends. However, as mentioned in §7.2, these trends are 
statistically not very significant. 

In Section 7.1, we found smaller dust obscuration in 
AGN-dominated systems and wondered if it was an evo- 
lutionary effect. We revisit Figure [2TH . taking into ac- 
count the interaction classes and nuclear separations. In 
Figure \37\ ) and c, we do find trends moving from regions 
Rl thru R2 to R3 (see summary in Table \TS\) : (a) close 
pairs or singles (separation < 1 kpc) occupy 36/^1/63% 
o] the total in R1/R2/R3; (b) middle interaction stages 
(Illa/b) decrease from R1-R3, occupying 68/4-3/29% of 
the total; and (c) late interaction stages (IVa/b) increase 
from 26/39/54% of the total in R1/R2/R3. These trends 
with morphology suggest that the importance of dust ex- 
tinction generally peaks during the intermediate stages of 
merger evolution (Illb/IVab) before dust gets destroyed 
or blown away during the last phase of the merger. A 
similar trend was found by Rossa et al. (2007) in the 
(smaller) sample of galaxies of the Toomre sequence. 

Even in the absence of projection effects and misiden- 
tifications, there are theoretical grounds for significant 
scatter in the simple evolutionary picture outlined above. 
Numerical simulations have been used extensively to 
study the dynamical evolution of merging galaxies and 
the associated inflow of gas to the central regions (e.g., 
Barnes & Hernquist 1996; Mihos & Hernquist 1996; 
Springel et al. 2005; Cox et al. 2006; Naab et al. 2006). 
These simulations capture star formation and AGN fuel- 
ing via "sub-resolution physics" : phenomenological mod- 
els which tie the star formation and AGN accretion - 
and their subsequent feedback on the gas - to the phys- 
ical properties of the gas on ~100 pc scales. Because 
the physical scales for accretion onto AGN are orders of 
magnitude smaller, these simulations can only provide a 
broad brush picture of the evolution of activity in merg- 
ing galaxies. Nonetheless, even with these limitations 
they provide a plausible framework for discussing both 
the trends and the scatter in this evolutionary picture. 

The relative strengths of the first and second inflow 
phases depend on a wide variety of factors, which will 
lead to significant scatter in the evolutionary paths of 
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mergers. As the early inflow is moderated by dynamical 
instabilities in the host galaxies' disks, this phase is very 
sensitive to the intrinsic properties of the host, such as 
the presence or absence of central bulges to stabilize the 
disks (Mihos & Hernquist 1996), the disk surface mass 
density (Mihos et al. 1997), or the gas fraction of the 
disk (Springel et al. 2005). Galaxies more susceptible to 
these disk instabilities are more likely to suffer early in- 
flow and onset of AGN activity, increasing the scatter in 
AGN properties in the early interaction stage. Other fac- 
tors playing into the merger evolution include the orbital 
geometry (e.g., prograde versus retrograde encounters) 
and mass ratio of the encounter. However, simulations 
show that these factors are secondary to the structural 
properties of the galaxies, as long as we consider ma- 
jor mergers like ULIRGs (e.g., Mihos & Hernquist 1996; 
Springel et al. 2005; Younger et al. 2008). 

The strengths of the starburst and AGN activity dur- 
ing different stages may also couple via the induced ac- 
tivity. The early fueling of starbursts and AGN can both 
deplete and heat the gas, potentially limiting the ability 
to form powerful starbursts or AGN late in the merging 
process. Simulations which include both AGN and star- 
burst heating (e.g., Springel et al. 2005) suggest this is 
a small effect early in the encounter when there is ample 
fuel supply, but that at later stages the AGN heating is 
sufficient to cease further starburst activity (DiMatteo 
et al. 2005). The fact that we see examples of post- 
mergers with significant starburst activity is somewhat 
problematic for these models, and may indicate that the 
AGN feedback models may be overly efficient in these 
simulations. Weaker feedback could halt accretion near 
the AGN without terminating star formation in the more 
extended distribution of gas. 

The fact that there are multiple inflow epochs along the 
merging sequence implies that a scenario in which AGN 
only turn on at the final stages of coalescence is over- 
simplified. Indeed, depending on the complex interplay 
of these factors, multiple bursts of strong starburst or 
AGN activity throughout the process are not completely 
unexpected. In a probabilistic sense, the likelihood of 
strong AGN fueling is highest as the galaxy nuclei coa- 
lesce, due to the rapidly varying gravitational potential 
that drives high inflow rates. However, the potential for 
AGN fueling exists throughout the interaction process. 
In a large sample like ours, this trend can be seen in the 
data, where later interaction classes and single nucleus 
objects do indeed have higher AGN fractions and larger 
Eddington ratios on average. However, the large scat- 
ter also indicates that AGN activity can in some cases 
dominate the galaxy's radiative energy output even at 
intermediate merger stages. 

Numerical simulations plausibly show that the evolu- 
tion of starburst and AGN activity among mergers can 
vary signiflcantly from merger to merger, due to varia- 
tions in the global properties of the encounter and the 
progenitor galaxies. However, there may be scatter in 
evolutionary paths that is unresolved by the simulations. 
Ultimately (and unfortunately) the detailed predictions 
for star formation and AGN accretion in the models de- 
pend critically on the sub-resolution physics. Varying the 
hydrodynamical equation of state or the prescriptions for 
star formation or AGN activity can significantly change 
the detailed results (e.g., Barnes 2004; Cox et al. 2006). 



Furthermore, the sub-resolution density structure of the 
gas is critical for driving continued inflow from the '--^lOO 
pc resolution limit of the models down to the accretion 
scale of the central AGN. Given these uncertainties and 
limitations inherent to the simulations, it is likely that 
the simulations merely give a time averaged expectation 
for AGN and starburst activity, averaged over the dy- 
namical timescale for the inner few hundred pc. The 
instantaneous rate of activity (as measured by our ob- 
servational data set) may show significant time variation 
- not seen in simulations - due to stochastic, small scale 
physical processes. This may also account for a signifi- 
cant amount of the scatter seen in the evolutionary trends 
shown in Figures 36 and 38. This random component of 
accretion is known to be important among local AGN 
(e.g., Davies et al. 2007). 

In summary, we see trends in AGN fraction, nuclear 
obscuration, and possibly Eddington ratio as a function 
of interaction class and nuclear separation, but with con- 
siderable scatter along the merger sequence. As shown in 
Figure \37\ ) and c, roughly half of fully merged ULIRGs 
have not (yet) succeeded in producing AGN-dominated 
systems and some pre-merger ULIRGs are already AGN- 
dominated. Some part of this scatter may be attributed to 
projection effects and misidentifications (that is, obser- 
vational effects). However, an equal or greater portion of 
the scatter is probably due to the physics identified by nu- 
merical simulations: the varying initial conditions among 
interactions and the fact that starburst and AGN activ- 
ity can peak locally in intensity prior to final coalesence. 
Further scatter may arise due to small-scale stochastic 
processes that are presently unresolved by these simula- 
tions. 

A revision of the evolution scenario of Sanders et al. 
(1988a, 1988b) is needed to explain all of these results. 
The "softer" version we propose requires the presence of 
multiple evolutionary paths in the phase space of AGN 
contribution, Eddington ratio, and dust obscuration ver- 
sus merger phase. For AGN contribution and Eddington 
ratio, this almost certainly includes paths that are not 
monotonically increasing with time. 

8. SUMMARY 

We have carried out a detailed Spitzer IRS study of 
the MIR continuum, absorption, and emission line prop- 
erties of a carefully selected sample of 74 ULIRGs and 34 
PG QSOs within z - 0.3. For the first time in ULIRGs, 
the continuum and dust features were modeled using a 
combination of PAH templates, blackbodies punctuated 
by deep extinction and absorption features, and silicate 
emission features, when necessary. The main observa- 
tional results are the foUowings: 

1. We find that the /30//15 and (PAH-free) MIR/FIR 
flux ratios are powerful continuum diagnostics of 
AGN activity among ULIRGs and QSOs. 

2. We confirm the broad range of silicate obscura- 
tion among ULIRGs, with the optically classified 
Seyfert ULIRGs being less obscured on average 
than the HH-like and LINER ULIRGs. The loose 
correlations seen between silicate, H2O ice -I- hy- 
drocarbons, C2H2, and HCN absorption features 
imply significant variations in composition of the 
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dense absorbing material from one ULIRG to the 
next. The HCO+ 12.1 fim and HNC 21.7 
features were not detected in any individual ob- 
ject or the average spectrum of the more obscured 
ULIRGs. 

3. The average PAH-to-FIR flux ratio of ULIRGs is 
remarkably similar to that of PG QSOs. No ob- 
vious trend is seen with optical spectral type, but 
both L(PAH)/L(FIR) and the 7.7/xm PAH equiv- 
alent width decrease with increasing extinction in 
HII/LINER ULIRGs. We confirm the strong cor- 
relation of EW(7.7/xm PAH) with optical spectral 
typos and find a similarly strong correlation with 
the /30/./15 ratio. 

4. Our analysis of the fine structure lines in ULIRGs 
and QSOs reveals a continuous excitation sequence 
with the cool {f25ixia/ feonm. < 0.1) optically clas- 
sified Hll-like and LINER ULIRGs at the low- 
excitation end of the sequence, the PG QSOs at the 
high-excitation end, and the warm optically classi- 
fied Seyfert ULIRGs in between. 

5. Warm H2 masses range from ^ 0.5 to 20 x 10^ 
Mq with an average (median) of ~ 3.8 (3.3) and 
3.6 (3.2) X 10^ M0 for the ULIRGs and QSOs, 
respectively. These masses arc typically a few per- 
cent of the cold gas mass derived from ^^CO obser- 
vations. The temperature-sensitive H2 S(2)/S(l), 
S(3)/S(l), and S(3)/S(2) flux ratios suggest pos- 
sible heating by the AGN in Seyfert ULIRGs and 
PG QSOs and shock excitation in LINER ULIRGs. 
However, dust extinction and/or variations in the 
ortho-to-para ratio make the results inconclusive. 

6. The average MIR spectnun of PAH-dominated 
ULIRGs suggests supersolar neon abundance while 
optical spectra indicate roughly solar oxygen abun- 
dance. Uncertainties on the exact value of the so- 
lar Ne/0 ratio may (partly) erase this discrepancy. 
However, if confirmed, this discrepancy may im- 
ply that these two methods trace different layers 
or star-forming regions of the ULIRGs, different in 
both extinction and abundance. No trend is seen 
in the sample galaxies between optically-derived 
metallicity and emission-line, absorption-line, and 
continuum properties. This result is not surprising 
given the relatively narrow range of optical metal- 
licity (0.6 - 2.4 X solar) covered by our sample and 
the small number (13) of objects for which we have 
reliable metallicity measurements. 

7. The contribution of an AGN to the bolomctric lu- 
minosity in these systems is quantified using six 
different methods based on (1) the [O IV] 25.9 
/im/[Ne II] 12.8 /im ratio, (2) the [Ne V] 14.3 
/zm/[Ne II] 12.8 /im ratio, (3) the equivalent width 
of PAH 7.7 fim, (4) the PAH (5.9 - 6.8 fj,m) to con- 
tinuum (5.1 - 6.8 /im) flux ratio combined with the 
continuum (14 15 /im) / (5.1 - 5.8 /im) flux ratio, 
(5) the MIR blackbody to FIR flux ratio, and (6) 
the /30//15 continuum flux ratio. Good agreement 
to within ~ ±10—15% on average is seen amongst 
the various methods. This agreement rules out the 



possibility that a MIR-buried but FIR-bright AGN 
is present in many of these objects. 

From these results we draw the following three main con- 
clusions: 

1. The average AGN contribution in ULIRGs is ^ 
35—40%, in agreement with previous ISO stud- 
ies. Strong correlations exist between AGN con- 
tributions, optical spectral types, and /25//6O ra- 
tios. The AGN contributions range from ~15— 35% 
among cool HII and LINER ULIRGs to r^50% and 
75% among warm Seyfert 2 and Seyfert 1 ULIRGs, 
respectively. The PG QSOs fall along the extrapo- 
lation of these trends, with AGN contributions typ- 
ically larger than ^ 80%. The largest AGN con- 
tributions are also observed at the smallest nuclear 
separations and latest interaction classes. 

2. All ULIRGs in our sample fall in three dis- 
tinct AGN classes: (1) Objects with small extinc- 
tions and large PAH equivalent widths are highly 
starburst- dominated, (2) System,s with large extinc- 
tions and modest PAH equivalent widths have larger 
AGN contributions, but are still mostly starburst- 
dominated, (3) ULIRGs with both small extinc- 
tions and PAH equivalent widths are either AGN 
dominated or show a balance between starburst and, 
AGN. The AGN contributions in highly obscured, 
class 2 ULIRGs are necessarily more uncertain 
than in the other objects, and we cannot formally 
rule out the possibility that these objects represent 
a physically distinct type of ULIRGs. However, a 
weak trend is seen toward smaller nuclear separa- 
tions and later merger stages along the sequence (1) 
- (2) - (3). These results suggest that d,ust extinc- 
tion generally peaks during the intermediate stages 
of merger evolution, before the dust gets destroyed 
or blown away during the late-merger phase. 

3. A "softer" version of the standard ULIRG - QSO 

evolution scenario is needed to explain the scatter 
in trends of AGN contribution, Eddington ratio, 
and dust obscuration with merger stage. With our 
large sample size we are able to discern the average 
trends discussed above. However, roughly half of 
fully merged, ULIRGs have not (yet) succeeded in 
producing AGN-dominated systems or blown away 
their obscuring dust screen, and some pre-merger 
ULIRGs are already AGN-dom,inated. Our revised 
evolutionary picture permits m,ultiple paths that are 
not necessarily monotonic in quantities like AGN 
contribution and Eddington ratio. Such a scenario 
is consistent with numerical simulations of merger- 
induced starburst and A GN activity. These simula- 
tions show the highest inflow rates when the galax- 
ies coalesce, but allow for significant episodes of in- 
flow and nuclear activity throughout a major galaxy 
merger. The strength and timing of these episodes 
will vary depending on the initial conditions of the 
interaction, and quite possibly on stochastic pro- 
cesses presently unresolved by simulations. 

Finally, we point out that the continuum-based meth- 
ods used here to quantify the power source in these local 
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ULIRGs and QSOs arc ideally suited to the study of 
faint high-z systems. In particular, the calibration of the 
/30//15 method we derived from our local sample may be 
used in the future to quantify the AGN contribution in 
U/LIRGs at z ~ 1 - 1.5 using the MIPS 70/24 /xm flux 
ratio (e.g., Sajina et al. 2007) and at higher redshifts 
with Herschel. 
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APPENDIX A. SIX METHODS TO DERIVE AGN CONTRIBUTIONS 
The six methods used to calculate the contribution of the AGN to the bolometric luminosity of the ULIRG or 
PG QSO are described briefly in the notes to Table [TTl Here we first discuss the assumptions that apply to all of the 
methods, and then describe each method individually with their respective strengths and weaknesses. 

For each of these methods, we compare the observed quantities derived from our data with pure- AGN and pure- 
starburst zero points [for method #3, we also compare the data with a pure-PDR (photodissociation region) zero 
point]. The pure- AGN zero point is set to the average value of the FIR- undetected (unobscured) PG QSOs in our 
sample, while the pure-starburst zero point is set to the average value of the most starburst-like ULIRGs in our sample. 
The results of these comparisons provide AGN fractional contributions to the [Ne II], 8 fvca. continuum, 5.3 - 5.8 /.tm 
continuum, FIR, and 15-/im luminosities, respectively. Next, we apply correction factors to transform these various 
luminosities into bolometric luminosities. These correction factors are listed in Table [10] for starburst (average of 
Hll-like ULIRGs) and AGN (average of FIR- undetected PG QSOs). The method used to calculate these correction 
factors is described in the note to that table. For this and the other methods, we assumed L{hol) = 1.15 i(IR) for 
aU ULIRGs (Kim & Sanders 1998) and L{bol) = 7 L(5100 k)+L{m) for aU PG QSOs. The latter includes both 
contributions from the 'intrinsic' AGN luminosity (Paper II), as well as AGN and starburst luminosity reprocessed by 
dust. 

We make no attempts to correct our data for extinction. The impact of dust extinction on the various features 
(continuum, fine structure lines, PAH features) depends greatly on the distribution of the dust relative to the sources 
of emission (e.g., dust screen versus mixed distribution) and we have very little information on this issue (except for the 
fact that it is almost certainly lower for PAHs than for the continuum and that foreground dust screens fit the continuum 
better than mixed dust screens; Sections 5.2 and 6.3). So, rather than making ad hoc assumptions on the dust geometry 
and running the risk of producing unphysical results (e.g., absurdly high PAH or MIR luminosities), we did not apply 
any dust extinction correction to the measured quantities. A comparison of our (extincted) measurements with those 
of (unobscured) FIR-undetected PG QSOs may therefore underestimate the AGN contribution in our objects. 

Dust extinction is not an issue when comparing our data with the pure-starburst zero point since this latter is based 
on the observed values of the most starburst-like ULIRGs. Starburst galaxies of lower infrared luminosities (lower star 
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formation rates and/or extinction) were not used to set this zero point or for the starburst bofometric corrections. 
Indeed, as pointed out by a number of studies, the unusual conditions (higher density, more intense radiation field) 
in ULIRGs produce systematic shifts in the relations between FIR, [Ne II] and PAH emission in such a way that we 
cannot use the relations derived from normal starburst galaxies to quantify ULIRGs. The competition of the dust with 
the gas for absorption of the ionizing photons becomes increasingly more effective as the density of the star-forming 
regions increases, so the [Ne IIJ/FIR ratio is reduced in ULIRGs (e.g., Rigby & Rieke 2004; Dopita et al. 2006; Calzetti 
et al. 2007). The more intense radiation field in ULIRGs induces a greater ionization or dehydrogenation of the PAHs, 
so the PAH/FIR ratio is also reduced (e.g., Tielens et al. 1999; Helou et al. 2001). Pure-starburst zero points and 
bolometric corrections based on starburst galaxies of lower infrared luminosities would therefore underestimate the 
true contribution of star formation in ULIRGs.^ 

Method #1 ([O IV]/[Ne II] ratio) and Method #2 ([Ne V]/[Ne II] ratio). The large number of ULIRGs 
without firm [O IV] and/or [Ne V] detection makes the results based on the individual [O IV] /[Ne II] and [Ne V]/[Ne II] 
ratios subject to potentially large systematic uncertainties. We treated upper limits as detections in our analysis so 
the AGN contribution as derived with this method should be considered upper limits as well. To further constrain 
these numbers, we measured the line ratios from average spectra, produced by normalizing individual spectra in each 
category listed in Table [TT] (e.g. spectral type, /25//6O1 infrared luminosity, morphological classes) to the same [Ne II] 
or FIR flux. The resulting average AGN fractional contributions, while affected by systematic uncertainties arising 
during the averaging procedure, are roughly consistent with the numbers in Table fTTl The average spectra suggest that 
the actual average AGN contributions are generally not lower by more than about 5—10% from the limits themselves, 
although for a few categories (e.g., H II regions), the actual values may be lower by a larger factor (10—20%). The 
pure-AGN zero point of the [O IV]/[Ne II] method is set at log([0 IV]/[Ne II]) = 0.6, corresponding to the average 
value of the FIR- undetected PG QSOs in our sample, while the pure-starburst ratio [O IV] /[No II] is set at zero (as 
described in Table [9l the non-zero ratios actually measured in starbursts are negligible for the purposes of computing 
AGN contributions). The assumption here is that [Ne II] emission from a pure starburst lowers [O IV] /[Ne II] below 
the pure AGN value, as shown by the diagonal line in Figure This figure also shows the calibration based on 
the [Ne V]/[Ne II] ratio. The pure-AGN log([Ne V]/[Ne II]) = 0.10 and is based once again on the average value of 
FIR-undetected PG QSOs. The correction factors used to transform [Ne II] luminosities into bolometric luminosities 
are listed in Table [TOl Note that Armus et al. (2007) did not apply this last correction to their numbers so their AGN 
fractional contributions relate to the [Ne II] luminosities, not the bolometric luminosities, and are considerably smaller 
than the numbers presented here. 

Method #3 (PAH 7.7 fim equivalent width). To facilitate comparisons with most of the published ISO results 
(e.g., Genzel et al. 1998; Lutz et al. 1999; Rigopoulou et al. 1999; Tran et al. 2001), we used the PAH 7.7 ^m equivalent 
widths to quantify the role of AGN in the 1-Jy ULIRGs. Recall that our fits allow only a small range 0.13 dex) in 
PAH 6.2/7.7 and 11.3/7.7 ratios so the conclusions based on the 7.7 fim feature also apply to first order to the other 
PAH features (see description of Method ^4 below). The pure-starburst zero point of this method, log[EW(PAH 7.7 
/im)] = 0.75, is near the maximum value observed in our sample. The AGN will reduce this quantity by contributing to 
the continuum emission at this wavelength and destroying the PAH molecules (e.g., Voit 1992). We assume EW(PAH 
7.7 ^m) = for a pure AGN and PAH destruction due to AGN radiation that is proportional to the AGN fractional 
contribution. The results do not depend sensitively on this last assumption. The correction factors used to transform 
8 /im continuum luminosities into bolometric luminosities are listed in Table 1101 

Method #4 (modified Laurent et aL method). This method is inspired by Laurent et al. (2000), but uses the 
modifications of Armus et al. (2007) to avoid contamination by PAHs in the continuum fluxes (5.3 - 5.8 /xm instead of 
5.1 - 6.8 fim). Figure [Ml shows the results for our sample. The zero points for the pure starburst and PDR are from 
Armus et al. (2007) and the zero point for the pure AGN corresponds to the average value for the FIR-undetected 
PG QSOs to reduce possible starburst contributions to the continuum emission (Paper II). We have moved the Armus 
et al. pure PDR point to the right by 0.3 dex to encompass the group of points that would otherwise fall outside the 
mixing region. Note that there is also uncertainty in the "pure AGN" point, since this is an average spectrum and the 
PAH luminosity is only an upper limit. The correction factors to transform 5.3 - 5.8 fim continuum luminosities into 
bolometric luminosities are listed in Table [TUl 

Method #5 (PAH-free 5 — 25 /im to FIR continuum ratio). In Section 6.1.3, we showed that the PAH-free, 
silicate-free MIR(5 - 25 /i)-to-FIR ratio derived from our fits was an excellent probe of nuclear activity in the 1-Jy 
ULIRGs (FigurelKl)- We adopt log[L(MIR)/L(FIR)] = 0.35 and -1.25 as the zero points for pure AGN and starburst, 
respectively. The AGN contribution is calculated from a linear interpolation between these two extremes. The zero 
point for the pure AGN corresponds to the average MIR/FIR ratio of FIR-undetected PG QSOs. The zero point for 
the pure starbursts is calculated from the ten ULIRGs with the lowest MIR/FIR ratios. The correction factors to 
transform FIR luminosities into bolometric luminosities are listed in Table [TUl 

Method #6 (/30//15 continuum ratio). The application of method ^5 to a large number of ULIRGs is time 
consuming since it involves detailed template fitting of the MIR SED. A more straightforward method, method ^6, is 
based on the /30//15 continuum ratio, which was found to be more tightly correlated with the PAH-free, silicate-free 

^ Another concern about the use of published [Ne IIJ/FIR and PAH/FIR ratios of normal starburst galaxies is aperture effects: the 
[Ne II] and PAH features are derived from Spitzer spectra with entrance apertures that are much smaller than IRAS, from which the FIR 
fluxes are derived. This effect will underestimate the actual [Ne IIJ/FIR and PAH/FIR ratios of normal starburst galaxies. In contrast, 
these aperture effects do not affect the ratios of ULIRGs significantly because most of the FIR, [Ne II], and PAH emission is produced 
within the central kpc of these objects, so is contained well within the Spitzer apertures. 
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MIR/FIR ratio than any other MIR continuum ratio at our disposal (Figure [TB)) . This method is roughly equivalent 
to using the MIPS 70/24 /xm flux ratios to search for AGN activity in z ~ 1 U/LIRGs. Here we use /30//15 as a 
surrogate of the PAH- free, silicate-free MIR/FIR ratio and adopt log(/3o//i5) = and 1.35 as the zero points for pure 
AGN and starburst, respectively. The AGN contribution is calculated from a linear interpolation between these two 
extremes. The zero point for the pure AGN corresponds to the average /30//15 ratio of FIR-undetected PG QSOs, 
while the zero point for the pure starbursts is calculated from the ten ULIRGs with the largest /30//15 ratios. The 
correction factors to transform 15-/im luminosities into bolomctric luminosities are listed in Table [Till 
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TABLE 1 
Sample 



Galaxy 
(1) 



z 

(2) 



log(L(6oO/L0) 
(3) 



Type 
(4) 



IC 

(5) 



NS 
(6) 



Ref 
(7) 



ULIRGs 



F00091-0738 
F00188-0856 
F00397-1312 
F00456-2904:SW 

F001S2-2721 



0.118 
0.128 

0.262 
0.110 

0.129 



12.36 
12.43 

12.96 
12.29 
12.09 



HII 
L 
HII 
HII 

L 



Illb 

V 

V 
Ilia 

Illb 



2.31 

< 0.34 

< 0.61 
22.80 

7.,39 



References. — 1 = Sakamoto et al 1999; 2 = Scoville et al. 2000; 3 = Beswick et aL 2001; 4 = Kim et al. 2002; 5 = Veilleux et al. 2006; 6 = 
Veilleux et al 2009 

Note. — Col(l); Galaxy name. Coordinate-based names beginning with "F" are sources in the IRAS Faint Source Catalog. Col. (2): Redshift. 
Col.(3): Bolometric luminosity. For ULIRGs, we assume L{hoi) = 1.15L(IR). For PG QSOs, we assume L{bol) = 7L(5100 A) + L{IR) (Netzer et 
al. 2007). Col. (4): Optical spectral type, from Veilleux et al. (1995, 1999a) and Rupke et al. (2005a). Col. (5): Interaction class, from Veilleux et 
al. 2009, Veilleux et al. 2006, or Veilleux et al. 2002 (in order of preference). Col. (6): Nuclear separation, in kpc. Col. (7): Reference for nuclear 
separation. 



TABLE 2 
Observations 



Galaxy 
(1) 



FID 

(2) 



SL2 
(3) 



SLl 
(4) 



Exposure Time 
SH 
(5) 



LH 

(6) 



LL 
(7) 



ULIRGs 



F00091-0738 
F00188-0856 
F00397-1312 
F00456-2904:SW 
F00482-2721 



3187 
105 
105 
3187 
3187 



240 
240 
240 
240 
240 



240 
240 
240 
240 
240 



480 
720 
720 
960 
720 



960 
480 
480 
720 
1440 



Note. — Col.(l): Galaxy name. Col. (2): Spitzer proposal lD(s) under which data was taken. Col. (3-7): Exposure times for each IRS module, 
in seconds. The LL exposure time is listed only when these data were used in the fit. 



TABLE 3 
Emission Line Fluxes 



Galaxy [NeVI]7.65 H2 S(3) 9.66 [S IV] 10.51 H2 S(2) 12.28 Hua 12.37 [NeII]12.81 [NeV114.32 [NeIII]15.55 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



H2 



ULIRGs 



F00091-0738 
F00188-0856 
F00397-1312 
F00456-2904:SW 
F00482-2721 



<5.18E-22 
7.27E-22( 8) 
3.89E-22(18) 
6.62E-22( 9) 
4.73E-22(10) 



<3.76E-22 
<8.96E-23 

2.66E-22(27) 
<1.24E-22 
<1.90E-22 



2.55E-22(59) 
2.99E-22(27) 
2.39E-22(38) 
4.14E-22(18) 
2.90E-22(15) 



<3.15E-22 
<8.99E-23 
<1.43E-22 
<9.59E-23 
<1.34E-22 



1.88E-21( 6) 
4.27E-21( 1) 
3.78E-21( 1) 
5.88E-21( 0) 
2.79E-21( 1) 



<3.92E-22 
<1.66E-22 
<2.56E-22 
<1.19E-22 
<1.55E-22 



<4.45E-22 
5.30E-22(10) 
2.14E-21( 9) 
1.50E-21( 2) 
6.04E-22( 6) 



<7. 
7. 
7. 
9. 



Note. — Col.(l): Galaxy name. Col.(2- 
in parentheses. Upper limits are 3(T. 



10): Atomic fine structure and H2 rotational emission line fluxes, in W cm . Percent errors are given 
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TABLE 4 
Emission Line Fluxes 



Galaxy 
(1) 



[Fell] 17.94 
(2) 



[S III] 18.71 
(3) 



[Ne V]24.32 
(4) 



[O IV] 25.89 
(5) 



[Fe II]25.99 
(6) 



H2 S(0) 28.22 
(7) 



[S III33.48] 
(8) 



[Si II]34.81 
(9) 



ULIRGs 



F00091-0738 <7.71E-22 



F00188-0856 
F00397-1312 
F00456-2904:SW 
F00482-2721 



<4.94E-22 
<5.89E-22 
<5.82E-22 
<1.69E-22 



<1.03E-21 

<3.09E-22 
1.29E-21(13) 
3.40E-21( 4) 
1.05E-21( 7) 



<9.27E-22 
<4.68E-22 
<5.52E-22 
<5.75E-22 
<2.75E-22 



<7.71E-22 
<4.92E-22 
<5.87E-22 
<5.82E-22 
6.75E-22(10) 



<7.71E-22 
<4.94E-22 
<5.89E-22 
<5.82E-22 
<1.69E-22 



<1.13E-21 
<5.76E-22 
<7.58E-22 
<6.21E-22 
<3.87E-22 



Note. — Col.(l): Galaxy name. Col. (2— 9): Atomic fine structure and H2 rotational emission line fluxes, in W cm ^. Percent errors are given 
in parentheses. Upper limits are 3cr. 



TABLE 5 
Continuum Measurements 



Galaxy 


6 /um 


15 //m 


20 nm 


25 fiiu 


30 fini 


12 ^.m (IRAS) 


25 ixm (IRAS) 


log[L(MIR)/L(FIR)] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




ULIRGs 


F00091-0738 


5.2 


117.7 


178.2 


546.1 


1145.5 


42.0 


325.0 


-0.73 


F00188-0856 


12.6 


78.3 


125.0 


327.8 


690.8 


26.5 


203.2 


-0.79 


F00397-1312 


78.3 


126.9 


146.8 


422.7 


792.8 


42.0 


157.8 


-0.30 


F00456-2904:SW 


8.9 


39.7 


111.8 


324.2 


691.9 


23.3 


210.7 


-0.98 


F00482-2721 


2.4 


9.3 


42.0 


129.9 


300.5 


5.2 


73.6 


-1.1 



Note. — Col.(l): Galaxy name. Col. (2-6): Rcst-framc Spitzer flux densities, computed from the IRS spectra using a 3.3% bandpass, in 
mjy. Col. (7-8): Observed-frame Spitzer fiux densities, computed using step function approximations to the IRAS 12 and 25 ^m system response 
functions. The flux given is the average /,y under the step function. Col. (9): Extincted 5 — 25 /^m luminosity minus PAH+silicate emission (i.e., 
blackbody only) as a (logarithmic) fraction of the far-infrared (40 — 122 jum) luminosity. For the PG QSOs, this latter quantity is available only 
for the three average spectra divided by L(FIR). 



TABLE 6 

1-"IT Hi:sn;rs: A.i;soni' i lox Air,Asrnic.Mi:xTs 



Galaxy t^^^^^ log[We<i(H20+HC)] logiVKeqCCzHa)] log[We<i(HCN)] 

(1) (2^ (3^ (4^ (5) 

F00091-0738 9.62 0.32 -2.09 -2.22 

F00188-0856 4.06 0.23 -2.84 -2.90 

F00397-1312 5.82 -0.81 

F00456-2904:SW 2.90 < 0.42 

F00482-2721 3.86 -0.21 



Note. — Col.(l): Galaxy name. CoI.(2): Effective peak silicate optical depth, computed using the ratio of the total extincted flux to the total 
unextincted flux. Col. (3-5): Rest-frame equivalent widths (in microns) of the water ice -t- hydrocarbon feature at 5 — 7 /^m; the C2H2 13.7 jLtm 
absorption feature; and the HCN 14.0 /^m absorption feature, respectively. 



TABLE 7 
Fit Results: PAH Measurements 



Galaxy 


\og[Weq{6^lm PAH)] 


log[Weq{7nm PAH)] 


Iog[L(PAH)/L(IR)] 


log[L(PAH)/L(FIR)] 


(1) 


(2) 


(3) 


(4) 


(5) 


F00091-0738 


-0.78 


-0.55 


-2.23 


-1.98 


F00188-0856 


-1.17 


-0.63 


-2.18 


-1.97 


F00397-1312 


-0.59 


0.11 


-0.83 


-0.59 


F00456-2904:SW 


0.14 


0.55 


-1.57 


-1.36 


F00482-2721 


-0.36 


0.14 


-2.04 


-1.84 



Note. — Col.(l): Galaxy name. Col. (2-3): Logarithmic rest-frame equivalent widths (in microns) of the PAH 6.2 jttm and 7.7 {^.m. features. 
Col. (4-5): Logarithmic ratios of the total PAH luminosity to the total infrared and far-infrared luminosities. 
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TABLE 8 
H2 Properties 



Galaxy 
(1) 



re:.(4-3) 

(2) 



Te:r(5-4) 

(3) 



log[M(H2)/M0] 
(4) 



ULIRGs 



F00091-0738 
F00188-0856 
F00397-1312 
F00456-2904:SW 
F00482-2721 



> 284.36 
305.07(14) 
274.51(19) 
317.97( 8) 
295.44( 8) 



> 391.26 
427.87(15) 
353.93(18) 
351. 26( 8) 
354.25( 8) 



8.34(31) 
9.16(45) 
8.37(19) 
8.46(18) 



Note. — Col.(l): Galaxy name. Col. (2 — 'Xj: Molecular hydrogen excitation temperatures determined from the S(l)/S(2) and S(2)/S(3) fluxes, 
respectively. Percent errors are given in parentheses. Col. (4): Molecular hydrogen mass computed using the partition function from Herbst et al. 
(1996), the S(l) flux, and the average of the (4 — 3) and (5 — 4) excitation temperatures. Percent errors are given in parentheses. 



TABLE 9 

Zero-Point Values for Computing AGN Contribution 



Quantity 


AGN 


SB/H II 


PDR 


(1) 


(2) 


(3) 


(4) 


log([0 IV] /[Ne II]) 


0.60 


... a- 




log([Nc V]/[Nc II]) 


0.10 


... a- 




log[Weq(PAH 7.7iim)] 




0.75 




log[/(PAH 6.2^lm)/f{5.3 - 5.8/im)] 


-1.82 


-0.40 


0.62 


log[/(14 - 16Mm)//(5.3 - 5.8//m)] 


0.07 


1.54 


-0.30 


Iog[L(MIR)/L(FIR)] 


0.35 


-1.25 




log[/3o//i5] 


0.20 


1.35 





Note. — Col.(l): Physical quantity by which AGN contribution to the bolometric luminosity is computed. Col. (2 — 4): Values of this quantity 
for a "pure" AGN, starburst / H II region, or photo-dissociation region (PDR). The values for [O IV]/[Ne II] and [Ne V]/[Ne II] are averages 
among the 8 PG QSOs undetected in the FIR (Netzer et al. 2007). The Weq(PAH 7.7ixm) value is the maximum value observed in our sample. The 
/(PAH 6.2;im)//(5.3 — 5.8;Am) and /(14 — 16;/m)//(5.3 — 5.8/Ltm) values are taken from Armus et al. 2007 for the H II region and PDR vertices, 
and from our FIR-undetected PG QSO subsamplc for the AGN vertex. Finally, the L(MIR)/L(FIR) and /30//15 values are estimated from our 
data, using the average of the FIR-undetected PG QSOs as a pure AGN. 

^ Technically, the [O IV] and [Ne V] emission in normal starbursts is non-zero (Lutz et al. 1998; Abel & Satyapal 2008); e.g., log([0 IV]/[Ne II]) 
~ -1.9 for the ISO starbursts with detected [O IV] (Verma et al. 2003). However, for the purposes of computing AGN contribution we can safely 
assume it is negligible. 



TABLE 10 

Bolometric Corrections for Computing AGN Contribution 



Quantity 


AGN 


SB 


(1) 


(2) 


(3) 


log([Ne II]/L(bo/)) 


-4.66 


-3.71 


log[L(5.3 - 5.8Atrn)/L(6oO] 


-1.92 


-2.85 


log[LA(8^im)/L(6o;)] 


-1.93 


-2.39 


\og[Lr,(lb^lnl)/L{bol)\ 


-14.33 


-14.56 


Iog[L(FIR)/L(6oO) 


-1.05 


-0.29 



Note. — Col. (2 — 3): Values of the quantity in column 1 for a "pure" AGN or starburst. Starburst values are averages over H II ULIRGs, 
and AGN values are averages over the 8 FIR-undetcctcd PG QSOs from Netzer et al. 2007. The bolometric luminosity is computed according to 
L{bol) = 1.15 L{IR) for all ULIRGs and L{bol) = 7 L(5100 A)+L(IR) for all PG QSOs. 
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TABLE 11 
Binned AGN Contribution 



Bin 

(1) 


No. 

(2) 


[O IV]/[NcII] 
(3) 


[Nc V]/[Nc II] 
(4) 


Wen (PAH 7.7fiin) 
(5) 


Laurent 

(6) 


L(MIR)/L(FIR) 

(7) 


/30//15 
(9) 


All 


All ULIRGs 


72 - 74 


38.5(it22.5) 


42.4(±26.3) 


40.0(±23.2) 


37.4(±26.5) 


31.5(±25.9) 


43.3(±27.1) 


38.8(±21.1) 


Spectral Type 


HII 
LINER 
Sey2 
Seyl 


17-18 
31-32 

12 

9 


<26.7(±15.2) 
<30.8(±12.9) 
58.5(±22.9) 
63.7(±27.9) 


<28.8(±19.6) 
<31.1(±17.4) 
70.2(±20.5) 
72.2(±23.6) 


29.1(±23.1) 
36.9(±21.0) 
48.8(±18.8) 
65.2{±14.1) 


20.4(±16.8) 
35.8(±23.8) 
44.8{±18.0) 
75.8(±22.5) 


22.5(±19.0) 
20.4(±18.0) 
40.7(±21.3) 
76.0(±17.4) 


33.1(±24.9) 
32.2(±19.6) 
55.7(±24.9) 
83.3(±12.3) 


27.1(±15.0) 
31.2(±14.6) 
53.1(±16.2) 
72.7(±14.4) 


log(/25//60) 


X < -1.2 
-1.2 <X < -1.0 
-1.0 <X < -0.8 

-0.8 < X 


9-10 
29-30 

13 

21 


<29.4(±14.4) 
<27.4(±12.5) 
43.8(±22.8) 
54.5(±26.0) 


<40.4(±18.2) 
<27.1(±17.0) 
44.8(±24.6) 
62.7(±28.0) 


32.9(±21.6) 
30.4(±19.2) 
39.4(±21.6) 
57.4(±21.1) 


28.3(±17.2) 
28.2(±19.3) 
37.2(±26.6) 
55.0(±31.2) 


9.5(± 8.7) 
16.4(±12.3) 
30.1(±15.1) 
64.2(±18.6) 


25.0(±18.5) 
27.7(±17.8) 
43.6(±21.8) 
74.3(±15.5) 


27.3(±12.5) 
26.4(±12.1) 
39.8(±16.2) 
61.3(±18.8) 


log(L(/R)/L0) 


X < 12.2 
12.2 < X < 12.4 
12.4 < X < 12.6 

12.6 < X 


27 
23 - 24 
14 
8-9 


<37.2(±21.4) 
<31.1(±18.0) 
43.7(±24.3) 
55.4(±27.1) 


<37.7(±27.4) 
<39.5(±24.2) 
45.3(±24.8) 
61.0(±27.3) 


34.7(±25.5) 
35.7(±21.9) 
46.0(±18.4) 
57.9(±16.6) 


32.7(±28.7) 
30.3(±19.4) 
44.5(±24.9) 
58.9(±28.8) 


27.8(±24.3) 
26.9{±21.3) 
33.3(±25.9) 
51.7(±34.9) 


38.5(±22.9) 
40.2(±28.8) 
49.0(±27.4) 
57.3(±31.8) 


34.8(±21.9) 
34.1(±16.1) 
43.6(±20.5) 
55.9(±24.4) 


Number of Nuclei 


Binary Systems 
Single Systems 


31 
34-36 


<36.5(±19.9) 
43.3(±25.0) 


<36.3(±24.5) 
51.3(±26.7) 


35.3(it25.6) 
45.3(±21.0) 


31.2(±25.4) 
43.5(±27.7) 


23.2(±21.3) 
40.4(±28.5) 


30.7(±25.1) 
54.6(±25.6) 


32.2(±20.3) 
46.2(±20.9) 


Nuclear Separation (kpc) 


X <1 
1< X <6 
X >6 


28 
14 
17 


<42.1(±25.1) 
<34.9(±18.9) 
<36.4(±22.0) 


<47.2(±28.1) 
<31.1(±22.9) 
<40.1(±26.1) 


44.0{±22.2) 
30.3(±25.8) 
38.4(±26.3) 


43.8(±29.8) 
28.1(±25.5) 
33.1(±26.2) 


41.2(±30.2) 
19.6(±20.1) 
26.9(±21.8) 


51.7(±28.2) 
25.9(±25.3) 
36.7(±23.8) 


45.0(±23.1) 
28.3(±20.2) 
35.3(±20.5) 


Interaction Class 


Ilia 
Illb 
IVa 
IVb 
V 


9 
22 
7-8 
18-19 
9 


<35.5(±23.9) 
<36.9(±18.7) 
<31.5(±14.6) 
48.9(±30.2) 
41.5{±17.3) 


<42.7(±26.1) 
<33.6(±23.9) 
<37.1(±16.0) 
58.6(±30.7) 
47.7{±20.6) 


37.4(±25.1) 
34.4(±26.4) 
36.5(±18.4) 
46.4(±22.7) 
50.9(±18.7) 


37.7(±30.8) 
28.5(±23.1) 
26.5(±15.4) 
47.8{±29.4) 
49.5(±29.0) 


26.1(±17.8) 
22.0(±22.9) 
15.1(±13.1) 
48.1(±26.5) 
47.0(±31.5) 


38.9(±17.1) 
27.4(±27.3) 
33.1(±13.6) 
60.4(±26.3) 
61.2(±23.5) 


36.4(±18.9) 
30.5(±21.0) 
29.7(± 8.2) 
51.6(±22.1) 
49.6(±20.0) 



Note. — Bolometric corrections arc computed using the denominator (N) in each quantity above, according to the formula 
AGN%(L(6oO)/100 = L{bol)''3" /{L{bol)''3" + L(boiy'') = 1/{1 + [100/AGiV%(A') - 1] x {N / L(bol))agn / {N / L{bol))sb}. Col.(l): Range of 
quantity over which AGN contribution is computed. Col. (2): Number of galaxies in each bin. Col. (3): Average AGN contribution to the 
bolometric luminosity computed from the [O IV]/[Ne II] line ratio, with standard deviation listed in parentheses. Individual upper limits are 
included in the average calculation, and those categories where upper limits dominate the average arc labeled as upper limits. AGN%/100 = 
[Ne II]ag„/([Ne ll]starburst + [Ne II]ogn) = [O IV]/[Ne ll\obseTvedK[0 IV]/[Ne II]agn - [O IV]/[Nc ll\starburst)- In the latter expression, 
we assume constant line ratios for a pure AGN or starburst. Col. (4): Same as column 3, but for [Ne V]/[Ne II]. Col. (5): Average AGN con- 
tribution computed from the equivalent width of the PAH T.Tfim feature, Weg(PAH 7.7^m), with standard deviation listed in parentheses. 

AGN% = fl^"^ continuum/I/l'' + /s ""1 = 1 " ^J W "j^'' {V AR 7.7^m)/H/|q''(PAH ^.^^ln\). In this calculation, wc assume PAH destruction due 

to AGN radiation, such that /o6s(PAH 7.7/im) = /s[,(PAH 7.7/^m) x (1 - AGN%). Col. (6); Average AGN contribution computed from the 
Laurent et al. 2000 diagram, as modified by Armus ct al. 2007, with standard deviation listed in parentheses. For 3-componcnt, 2-ratio 
mixing between an AGN, H II region, and PDR using quantities A = /(PAH 6.2/ini)//(5.3 — 5.8^m) and B = /(14 — 16^m) //(5.3 — 5.8/im), 
AGN%/100 = /a9n(5.3-5.8Mm)/[/a9n(5.3-5.8Atm) + /h2(5.3-5.8Mm)+/p<ir(5.3-5.8Mm)] = (AobsBh2+ApdrBobs+Ah2Bpdr-Aob8Bpdr- 
Ah2Bobs~ApdrBh2)/{AagnBh2+ApdrBagn+Ah2Bpdr — AagnBpdr — Ah2Bagn — ApdrBh2)- Col.(7): Average AGN Contribution computed 
from the L(MIR)/L(FIR) luminosity ratio, with standard deviation listed in parentheses. The PAH and silicate emission have been removed 
from the measured MIR luminosity; only the blackbody dust emission remains. AGN%/100 = L{FI R)agn/[L(FI R)i,b + L{FI R)agn] = 
[L{MIR)/L(FIR)„bs - L{MIR)/L{FIR)^b]/[L{MIR)/L{FIR)agn - L{MIR)/L{FIR)^b]- Col.(8): Average AGN contribution com- 
puted from the /3o/im//i5^m flux density ratio, with standard deviation listed in parentheses. The formula used is the same as for 
the L(MIR)/I/(FIR) diagnostic. Col. (8) Average-of-averages, with standard deviation listed in parentheses. We first average over the 6 
methods for each galaxy, and then average over all galaxies. [See Appendix A for more information on the individual methods.] 
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TABLE 12 
AGN Contributions 



Galaxy 
(1) 



1 

(2) 



2 

(3) 



Method 



3 

(4) 



4 

(5) 



5 

(6) 



6 

(7) 



Avg 
(8) 



ULIRGs 



F00091-0738 

F00188-0856 
F00397-1312 
F00456-2904:SW 
F00482-2721 



< 50.6 

< 21.0 

< 26.6 

< 18.5 
36.6 



< 63.9 

< 22.1 

< 33.6 

< 12.7 

< 29.2 



54.6 
57.4 
27.6 
8.1 
26.1 



16.6 
73.1 
44.1 
6.6 
29.4 



47.8 
52.5 
67.1 
15.6 
0.0 



26.1 
22.4 
59.0 
11.2 
5.0 



43.3 
41.4 
43.0 
12.1 
21.0 



Note. — Col.(l): Galaxy name. Col. (2 — 8): Percent of the bolometric luminosity produced by the AGN, as determined using six different 
methods, as well as the average percent over all six methods. Methods used: 1 = [O IV]/[Ne II]; 2 = [Ne V]/[Ne II]; 3 = IVeqCPAH 7.7(Um); 4 = 
Laurent; 5 = Z/(MIR)/L(FIR); and 6 = /so/Zis- [See Appendix A for more information on the individual methods.] 



TABLE 13 

Statistics on Regions of Weq(PAH 7.7/iM) vs. rg.y Space 







# of Galaxies 






% of Total 




Quantity 


Rl 


R2 




R3 


Rl 


R2 


R3 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 










Regions 


Subdivided by AGN % 









<20% 


15 


2 65% 


8% 


0% 


20 - 40% 


8 


12 1 35% 


48% 


4% 


40 - 60% 





11 11 0% 


44% 


44% 


>60% 





13 0% 


0% 


52% 


All 


23 


25 25 100% 


100% 


100% 


Regions Subdivided by Nuclear Separation 


>6 kpc 


8 


4 5 36% 


24% 


26% 


1 — 6 kpc 


6 


5 2 27% 


29% 


11% 


<1 kpc 


8 


8 12 36% 


47% 


63% 


All 


22 


17 19 100% 


100% 


100% 


Regions Subdivided by Interaction Class 


Illab 


13 


10 7 68% 


43% 


29% 


IVab 


5 


9 13 26% 


39% 


54% 


V 


1 


4 4 5% 


17% 


17% 


All 


19 


23 24 100% 


100% 


100% 



Note. — Col.(l): Quantity by which regions Rl - R3 in VKeg(PAH 7.7/Ani) vs. T9.7 phase space are subdivided. Col. (2 — 4): Number of galaxies 
in each region. Col. (5 — 7): Percentage of galaxies in each region. 
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TABLE 14 

Binned Eddington Ratio, Based on Photometric Black Hole Mass Estimates 



Bin 


No. 


[O IV]/[Ne II] 


[Ne V]/[Ne II] Weq(PAH 7.7fim) 


Laurent 


L(MIR)/L(FIR) 


/30//15 


All 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(9) 




All ULIRGs 


58 


-61 


-1.07(±0.38) 


-1.03(±0.41) 


-1.04(±0.47) 


-1.11 (±0.48) 


-1.21(±0.48) 


-0.97(±0.39) 


-1.07(±0.40) 


Spectral Type 


HII 


11 


- 12 


-1.09(±0.42) 


-1.03(±0.50) 


-1.01(±0.62) 


-1.15(±0.58) 


-1.11(±0.60) 


-0.83(±0.43) 


-1.03(±0.50) 


LINER 


24 


-26 


-1.06(±0.39) 


-1.07(±0.42) 


-1.00(±0.41) 


-1.09(±0.45) 


-1.32(±0.46) 


-1.02(±0.37) 


-1.10(±0.37) 


Sey2 




11 


-0.95(±0.30) 


-0.86(±0.29) 


-1.05(±0.32) 


-1.08(±0.31) 


-1.16(±0.30) 


-1.01(±0.36) 


-1.02(±0.27) 


Seyl 


8 


-9 


-1.12(±0.41) 


-1.04(±0.40) 


-0.97(±0.39) 


-0.91(±0.43) 


-1.00(±0.47) 


-0.95(±0.46) 


-1.03(±0.44) 


log(/25//60) 


X < -1.2 


7 


- 8 


-1.04(±0.39) 


-0.87(±0.34) 


-1.05(±0.65) 


-1.05(±0.47) 


-1.42(±0.42) 


-1.07(±0.39) 


-1.06(±0.38) 


-1.2 <X < -1.0 


21 


- 25 


-1.11(±0.36) 


-1.13(±0.42) 


-1.04(±0.50) 


-1.14(±0.50) 


-1.35(±0.49) 


-1.00(±0.37) 


-1.14(±0.42) 


-1.0 <X < -0.8 




11 


-1.03(±0.30) 


-1.00(±0.29) 


-1.06(±0.39) 


-1.12(±0.48) 


-1.22(±0.41) 


-1.05(±0.40) 


-1.08(±0.31) 


-0.8 < X 


16 


- 18 


-1.05(±0.46) 


-0.98(±0.48) 


-1.01(±0.42) 


-1.11(±0.50) 


-0.92(±0.42) 


-0.86(±0.41) 


-0.99(±0.44) 


\og{L{IR)/LQ) 


X < 12.2 


15 


- 17 


-1.10(±0.42) 


-1.10(±0.48) 


-1.16{d=0.40) 


-1.34(±0.55) 


-1.28{±0.42) 


-1.08(±0.37) 


-1.16(±0.42) 


12.2 < X < 12.4 


20 


- 23 


-1.12(±0.35) 


-1.05(±0.40) 


-1.10(±0.58) 


-1.14(±0.49) 


-1.22{±0.54) 


-0.92(±0.42) 


-1.09(±0.12) 


12.4 < X < 12.6 




13 


-0.98(±0.30) 


-0.96(±0.28) 


-0.94(±0.35) 


-0.98(±0.32) 


-1.17(±0.41) 


-0.92(±0.29) 


-0.99(±0.28) 


12.6 < X 


8 


-9 


-0.98(±0.51) 


-0.94(±0.51) 


-0.79(±0.36) 


-0.83(±0.36) 


-1.10(±0.57) 


-0.98(±0.50) 


-0.98(±0.47) 


Number of Nuclei 


Binary Systems 


21 


- 25 


-1.09(±0.34) 


-l.ll(it0.39) 


-1.14(±0.45) 


-1.24(it0.44) 


-1.40(±0.49) 


-1.11(±0.39) 


-1.18(±0.39) 


Single Systems 


34 


-35 


-1.04(±0.41) 


-0.96(±0.43) 


-0.96(±0.49) 


-1.03(±0.51) 


-1.07(±0.44) 


-0.89(±0.38) 


-0.99(±0.40) 


Nuclear Separation (kpc) 


X <1 


21 


-23 


-0.99(±0.41) 


-0.92(±0.45) 


-0.91(±0.44) 


-0.99(±0.55) 


-0.98(±0.45) 


-0.85(±0.40) 


-0.95(±0.42) 


KX <6 


10 


- 12 


-0.95(±0.29) 


-1.05(±0.37) 


-1.02(±0.38) 


-1.20(±0.49) 


-1.35(±0.57) 


-1.09(±0.51) 


-1.11(±0.41) 


X >6 


11 


- 13 


-1.21(±0.34) 


-1.17(±0.42) 


-1.22(±0.50) 


-1.27(±0.40) 


-1.44(±0.42) 


-1.14(±0.28) 


-1.25(±0.37) 


Interaction Class 


Ilia 




8 


-1.18(±0.32) 


-1.08(±0.44) 


-1.14(±0.42) 


-1.17(±0.36) 


-1.34(±0.30) 


-1.13(±0.28) 


-1.17(±0.32) 


Illb 


13 


- 17 


-1.04(±0.35) 


-1.13(±0.38) 


-1.14(±0.48) 


-1.27(±0.48) 


-1.43(±0.56) 


-1.10(±0.46) 


-1.19(±0.42) 


IVa 


7 


-8 


-1.13(±0.46) 


-1.07(±0.46) 


-1.14(±0.67) 


-1.25(±0.59) 


-1.44(±0.46) 


-1.11(±0.36) 


-1.19(±0.42) 


IVb 


17 


- 18 


-1.11(±0.39) 


-1.03(±0.44) 


-1.01(±0.42) 


-1.04(±0.45) 


-1.06(±0.37) 


-0.93(±0.35) 


-1.02(±0.38) 


V 




9 


-0.81(±0.36) 


-0.76(±0.38) 


-0.72(±0.37) 


-0.80(±0.50) 


-0.81(±0.41) 


-0.64(±0.35) 


-0.76(±0.37) 



Note. - Col.(l): Bin in which Eddington ratio is averaged. Col. (2): Number of galaxies in each bin. Col. (3-9): Average Eddington 
ratios, computed from the measured AGN luminosity and the photometrically-determined black hole mass. The AGN luminosity is 
measured as a fraction of the bolometric luminosity; this AGN contribution is estimated from each of 6 different mid-infrared diagnostics. 
Col. (10): Eddington ratio averaged over the 6 different ways of computing the AGN contribution. 



27 



TABLE 15 

Binned Eddington Ratio, Based on Dynamical Black Hole Mass Estimates 



Bin 


No. 


[O IV]/[Ne II] 


[Ne V]/[Nc II] 


Weq(PAH 7.7Atm) 


Laurent 


L(MIR)/L(FIR) 


/30//15 


All 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(9) 




All ULIRGs 


22 


-25 


-0.34(±0.61) 


-0.33(±0.67) 


-0.31(±0.66) 


-0.42(±0.69) 


-0.33(±0.70) 


-0.22(±0.67) 


-0.35(±0.63) 


Spectral Type 


HII 


7 


- 9 


-0.34(±0.63) 


-0.36(±0.70) 


-0.25(±0.73) 


-0.43(±0.76) 


-0.18(±0.72) 


0.01(±0.67) 


-0.34(±0.67) 


LINER 


5 


-6 


-0.54(±0.66) 


-0.60(±0.76) 


-0.57(±0.66) 


-0.70(±0.59) 


-0.63(±0.77) 


-0.61(±0.64) 


-0.61(±0.63) 


Sey2 


4 


-5 


-0.34(±0.63) 


-0.23(±0.65) 


-0.33(±0.82) 


-0.32(±0.81) 


-0.58(±0.79) 


-0.40(±0.78) 


-0.38(±0.70) 


Seyl 




4 


0.02(±0.62) 


0.07(±0.58) 


0.07(±0.51) 


0.09(±0.50) 


0.13(±0.47) 


0.17(±0.46) 


0.09(±0.51) 


log(/25//60) 


X < -1.2 


1 


- 3 


-0.59(±0.88) 


-0.46(±0.96) 


-0.25(±1.16) 


-0.30(±1.16) 


-0.00(±0.00) 


0.45(±0.00) 


-0.51(±0.83) 


-1.2 <X < -1.0 


5 


- 6 


-0.29(±0.68) 


-0.37(±0.77) 


-0.34(±0.74) 


-0.45(±0.81) 


-0.34(±0.96) 


-0.32(±0.89) 


-0.36(±0.78) 


-1.0 <X < -0.8 


4 


- 5 


-0.48(±0.50) 


-0.53(±0.53) 


-0.48(±0.67) 


-0.66(±0.50) 


-0.71(±0.70) 


-0.50(±0.68) 


-0.55(±0.54) 


-0.8 < X 


10 


- 11 


-0.23(±0.60) 


-0.19(±0.66) 


-0.23(±0.63) 


-0.33(±0.69) 


-0.17(±0.60) 


-0.11(±0.56) 


-0.22(±0.60) 


log{L{IR)/LQ) 


X < 12.2 


11 


- 13 


-0.55(±0.48) 


-0.60(±0.54) 


-0.58{±0.56) 


-0.69(±0.60) 


-0.65{±0.64) 


-0.54{±0.58) 


-0.60(±0.52) 


12.2 < X < 12.4 


5 


- 6 


-0.37(±0.53) 


-0.24(±0.51) 


-0.29(±0.54) 


-0.47(±0.41) 


-0.20(±0.50) 


-0.06(±0.48) 


-0.33(±0.46) 


12.4 < X < 12.6 


1 


- 2 


-0.25(±1.32) 


-0.37(±1.46) 


-0.18(±1.26) 


-0.21(±1.29) 


0.65(±0.00) 


0.73(±0.00) 


-0.25(±1.33) 


12.6 < X 




4 


0.36(±0.37) 


0.42(±0.41) 


0.36(±0.46) 


0.27(±0.64) 


0.26(±0.62) 


0.37(±0.64) 


0.34(±0.51) 


Number of Nuclei 


Binary Systems 


4 


- 6 


-0.36(±0.30) 


-0.48(it0.43) 


-0.43(±0.53) 


-0.70(it0.35) 


-0.53(±0.29) 


-0.53(±0.50) 


-0.50(±0.34) 


Single Systems 


14 


- 15 


-0.16(±0.64) 


-0.07(±0.61) 


-0.10(±0.63) 


-0.19(±0.69) 


-0.13(±0.77) 


-0.00(±0.69) 


-0.14(±0.65) 


Nuclear Separation (kpc) 


X <1 


16 


- 19 


-0.32(±0.68) 


-0.29(±0.74) 


-0.28(±0.70) 


-0.36(±0.74) 


-0.24(±0.80) 


-0.18(±0.76) 


-0.32(±0.71) 


KX <6 




1 


0.15(±0.00) 


0.25(±0.00) 


0.18(±0.00) 


-0.33(±0.00) 


-0.14(±0.00) 


0.13(±0.00) 


0.04(±0.00) 


X >6 


3 


-5 


-0.48(±0.18) 


-0.61(±0.24) 


-0.65(±0.41) 


-0.83(±0.31) 


-0.62(±0.20) 


-0.43(±0.30) 


-0.57(±0.21) 


Interaction Class 


Ilia 


1 


-2 


-0.56(±0.23) 


-0.80(±0.13) 


-1.08(±0.00) 




-0.69(±0.36) 


-0.47(±0.47) 


-0.68(±0.37) 


Illb 


3 


-4 


-0.26(±0.30) 


-0.33(±0.45) 


-0.21(±0.37) 


-0.55(±0.20) 


-0.43(±0.26) 


-0.56(±0.58) 


-0.41(±0.34) 


IVa 




3 


-0.43(±0.60) 


-0.45(±0.58) 


-0.40(±0.57) 


-0.76(±0.51) 


-0.57(±0.63) 


-0.38(±0.39) 


-0.50(±0.51) 


IVb 


7 


-8 


-0.33(±0.66) 


-0.19(±0.60) 


-0.31(±0.62) 


-0.31(±0.62) 


-0.28(±0.82) 


-0.19(±0.74) 


-0.31(±0.64) 


V 




4 


0.38{±0.35) 


0.47(±0.35) 


0.47(±0.31) 


0.44(±0.49) 


0.45(±0.52) 


0.61(±0.43) 


0.47(±0.39) 



Note. — Col.(l): Bin in which Eddington ratio is averaged. Col. (2): Number of galaxies in each bin. Col. (3-9): Average Eddington 
ratios, computed from the measured AGN luminosity and the dynamically determined black hole mass. The AGN luminosity is measured as 
a contribution of the bolometric luminosity; this AGN contribution is estimated from each of 6 different mid-infrared diagnostics. Col. (10); 
Eddington ratio averaged over the 6 different ways of computing the AGN contribution. 



TABLE 16 

Eddington Ratios from Photometry and Dynamics 



log{LAGN/LEdd), photometry \og{LAGN / LEdd), dynamics 

Galaxy 1 2 3 4 5 6 Avg 1 2 3 4 5 6 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



ULIRGs 



F00091-0738 


< 


-0.80 


< 


-0.70 


-0.77 


-1.29 


-0.83 


-1.09 


-0.91 


< 0.15 


< 0.25 


0.18 


-0.33 


0.13 


-0.14 


F00188-0856 


< 


-1.19 


< 


-1.17 


-0.75 


-0.65 


-0.79 


-1.16 


-0.95 














F00397-1312 


< 


-0.63 


< 


-0.53 


-0.61 


-0.41 


-0.23 


-0.29 


-0.45 


< 0.85 


< 0.95 


0.86 


1.07 


1.25 


1.19 


F00456-2904:SW 


< 


-1.41 


< 


-1.57 


-1.76 


-1.85 


-1.48 


-1.62 


-1.62 


< -0.73 


< -0.89 


-1.08 


-1.17 


-0.80 


-0.94 


F00482-2721 




-0.89 


< 


-0.99 


-1.04 


-0.98 




-1.75 


-1.13 















Note. — Col.(l): Galaxy name. Col. (2 — 8): log Eddington ratio, using AGN luminosity computed from six different methods and black hole 
mass from galaxy photometry, as well as the average ratio over all six methods of computing the AGN luminosity. Methods used: 1 — [Q IV]/[Ne II]; 
2 ^ [Ne V]/[Ne II]; 3 ^ W^^{PAU 7.7/xm); 4 ^ Laurent; 5 ^ L(MIR)/L(FIR); and 6 ^ /3o//i5- Col. (9 - 15): log Eddington ratio, using AGN 
luminosity computed, from the same six methods and black hole mass from galaxy dynamics, as well as the average ratio over all six methods of 
computing the AGN luminosity. 
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Fig. 1. — Distributions of (a) redshifts, (fe) infrared luminosities, and (c) 25-to-60/^m flux ratios of the Spitzer 1 Jy ULIRGs (black 
hatched histogram) vs. entire 1 Jy sample (blue histogram). In panels (b) and (c), IRAS 12 and 25/xm fluxes are used for the entire 1 
Jy sample. Many of these fluxes are upper limits, which we label explicitly in panel (c). IRAS-type Spitzer 12 and 25^m fluxes are used 
for the Spitzer subsample. The new flux measurements account for the imperfect overlap in panels (fe) and (c). The Spitzer ULIRGs are 
representative of the 1-Jy sample in both range and distribution of properties. 
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Fig. 2. — Distributions of (a) redshifts, (b) B-band absolute magnitudes, and (c) 25-to-60 ^tm flux ratios of the Spitzer PC QSOs (black 
hatched histogram) vs. entire PG QSO sample (blue histogram). As in Figure [1] IRAS 25/im fluxes are used for the entire sample, and 
IRAS-type Spitzer fluxes for the current subsample. The Spitzer QSOs sample the low redshift and low luminosity ends of the PG QSO 
sample. 
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0.0 0.1 0.2 0.3 0.4 



z 

Fig. 3. — Redshift distributions of the Spitzer VhlRGs (blue histogram) and PG QSOs (black hatched histogram). The two samples are 
well matched in redshift. 




log[ f,2 (IRAS) / Jy ] log[ fj^ (IRAS) / Jy ] 



Fig. 4. — Comparisons of IRAS (a) 12 and (b) 25 fim flux densities with Spitzer-derived quantities. Squares represent ULIRGs and 
circles are QSOs. Open symbols are upper limits. The colors of the symbols reflect the optical spectral types: red, green, blue, and black 
squares are Hll-like, LINER, Seyfert 2, and Seyfert 1 ULIRGs, respectively. The solid line is perfect agreement, and the dotted line is for 
f {I RAS)/ f (Spitzer) = 2. Excellent agreement is seen at the 5% level at 25 fim and at the 25% level at 12 /im. 
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Fig. 5. — Four examples of spectral decompositions: (a) PAH-dominated ULIRG, F15206+3342, (6) Highly obscured ULIRG, 
F00091-0738, (c) AGN-dominated ULIRG without siUcate emission, F11119+3257, and (d) AGN-dominated ULIRG/PG QSO with 
silicate emission, Mrk 1014 = PG 0157+001. In each case, the top panel shows the fit to the IRS spectrum, while the bottom panel shows 
the data-to-model flux ratio. The data are in black, the overall fit is in red, the three blackbody components are in green, the two PAH 
templates are in blue, and the silicate emission component is the long-dash blue line. See Section 5.2 for more detail on the fits. 
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Fig. 6. — IRS spectra of all ULIRGs considered in our study, in order of increasing right ascension. The sources of the overlaid photometry, 
shown as red stars, are discussed in Section 6. F11223— 1244:W is shown in this figure but was mistakenly omitted from the analysis. 
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5 6 7 8 9 10 20 30 

Wovelength (jim) 

Fig. 7. — Average IRS spectra for ULIRGs with 25-to-60 fim flux ratios, /25//6O1 above and below 0.1, compared with the QSOs in 
our sample (Paper II). The individual spectra in each category were normalized to have the same rest-frame 15 iim flux density. Note the 
progression from cool ULIRGs to warm ULIRGs, and then to QSOs. 
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5 6 7 8 9 10 20 30 

Wovelength (/xm) 

Fig. 8. — Average IRS spectra for ULIRGs with infrared luminosities larger or smaller than lO^'^ '^ Lq, compared with the QSOs in our 
sample. The individual spectra in each category were normalized to have the same rest-frame 15 fim flux density. 
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5 6 7 8 9 10 20 30 

Wovelength (/xm) 

Fig. 9. — Average IRS spectra for ULIRGs of various optical spectral types, compared with the QSOs in our sample. The individual 
spectra in each category were normalized to have the same rest-frame 15 fim flux density. Note the similarity between the average spectrum 
of Seyfert 1 ULIRGs and that of QSOs. 
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5 6 7 8 9 10 20 30 

Wovelength (/xm) 

Fig. 10. — Average IRS spectra for ULIRGs with PAH 7.7 fim equivalent widths larger or smaller than 1 fj,m, compared with the QSOs 
in our sample. The individual spectra in each category were normalized to have the same rest-frame 15 ^m flux density. 
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Fig. 11. — Average IRS spectra for ULIRGs with effective optical depth of the 9.7 fim absorption feature larger or smaller than 3.68 
(the sample median), compared with the QSOs in our sample. The individual spectra in each category were normalized to have the same 
rest-frame 15 ^im flux density. Significant PAH emission is detected in both average ULIRG spectra. 
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Fig. 12. — MIR color-color diagrams: (a) 15-to-6 vs. 30-to-6 fira flux ratios, (fe) 15-to-6 vs. 30-to-15 /xm flux ratios, (c) 25-to-60 vs. 
30-to-15 /im flux ratios. The meaning of the ULIRG and QSO symbols is the same as in Figure |3] In addition, the red stars and black 
triangles are starburst and Seyfert galaxies observed with ISO (Verma et al. 2003; Sturm et al. 2002; Brandl et al. 2006). The tightest 
correlation is seen in (c). 
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Fig. 13. — 5piteer-derived MIR (5 - 25 ^m) luminosities of 1-Jy ULIRGs and PG QSOs vs. FIR luminosities. The meaning of the 
symbols is the same as in Figure |4] The solid line is the line of equality, while the dashed lines show the locations of objects with FIR 
luminosities equal to (1/4, 1/2, 2, 4, 8, 16) X the MIR luminosities. All ULIRGs, except most of those that are optically classified as Seyfert 
Is, are MIR underluminous relative to QSOs. Seyfert 2 ULIRGs are intermediate between QSOs/Seyfert 1 ULIRGs and HII-like/LINER 
ULIRGs. The Hll-like ULIRG at very high infrared luminosity is F00397— 1312, and has heavily extincted PAH emission. 
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Fig. 14. — Distributions of the blackbody temperatures for the (a) cold, (b) warm, and (c) hot components of the template fits to the IRS 
spectra of the 1-Jy ULIRGs. The arrowhead in panel (c) indicates that the somewhat uncertain ~1000 K temperature are possibly lower 
limits (with an estimated range ~ 700 — 2000 K). The temperatures of the hot components in Hll-like / LINER ULIRGs are distinctly 
lower than in Seyfert ULIRGs. 




Fig. 15. — Distributions of the ratios of the monochromatic blackbody (i.e., excluding PAH emission) (a) 8, (b) 15, an d (c) 30 fim 
luminosities to the total infrared luminosities for 1-Jy ULIRGs of various optical spectral types (see the legend in Figure [T4j. The ratios 
involving the 8 and 15 fj,m luminosities are distinctly larger among Seyfert ULIRGs than among Hll-like / LINER ULIRGs. 



39 




1.4 : 
1.2 : 

^ 1.0 : 

to 

^ 0.8 : 

f 0-6: 
0.4 : 

0.2 : 
0.0 ^ 

1.5 - 



1'" 



en 
o 



(b) 



0.5 



0.0 I 



• .11 




(d) 



'I 



-1.5 -1.0 -0.5 0.0 0.5 -1.5 -1.0 -0.5 0.0 0.5 

log[ L(BB. 5-25/im, observed) / L(FIR) ] 

Fig. 16. — (a) Distributions of the ratios of 5 - 25 fim blackbody (i.e., ex clud ing PAH emission) luminosities to the FIR luminosities 
for 1-Jy ULIRGs of various optical spectral types (see the legend in Figure [M}. This ratio is distinctly larger among Seyfert ULIRGs 
than among Hll-like / LINER ULIRGs. (6, c, d) Ratios of 5 - 25 ^tm blackbody (i.e., excluding PAH emission) luminosities to the FIR 
luminosities for 1-Jy ULIRGs of various optical spectral types vs. 15-to-6 ^m flux ratios (6), 30-to-6 /im flux ratios (c), and 30-to-15 /im 
flux ratios (d). The meaning of the symbols is the same as Figure!?] addition, the small, medium-size, and large black circles correspond 
to the FIR-undetected, FIR-faint, and FIR-bright PG QSOs, respectively, as defined in Paper II. The strongest correlation is seen in (d). 
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Fig. 17 . — Distributions of the eflfective 9.7 /im silicate optical depth for 1-Jy ULIRGs of various optical spectral types (see the legend 
in Figure [T4t . The four Seyfert 1 ULIRGs with silicate in emission are not shown in this figure. Optically classified Seyferts generally have 
smaller optical depths than HII-like/LINER ULIRGs. 
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Fig. 18. — Effective silicate optical depths vs. (a) ratios of monochromatic blackbody 8 fim luminosities to total infrared luminosities, 
(fe) ratios of monochromatic blackbody 15 iim luminosities to total infrared luminosities, (c) 30-to-6 ^im flux ratios, and (d) 30-to-15 fim 
flux ratios for the 1-Jy ULIRGs of various optical spectral types. The meaning of the symbols is the same as in Figure |4] Seyfert ULIRGs 
generally populate the lower left portion of these diagrams. 
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Fig. 19. — Comparisons of absorption features: Effective optical depths of 9.7 /im silicate, equivalent widths of 14 fim HCN; equivalent 
widths of H2O ice + hydrocarbons; equivalent widths of 13.7 fim C2H2, and equivalent widths of 14 fim HCN. The meaning of the symbols 
is the same as in Figure |4] The strengths of these features are only loosely correlated, implying significant variations in the composition of 
the dense absorbing material from one ULIRG to the next. 



43 




Fig. 20. — {a),(c) Distributions of th e ra tio of the total PAH luminosities to FIR luminosities of the 1-Jy ULIRGS according to optical 
spectral types (see the legend in Figure I14I I and effective silicate optical depth (relative to the median value) . There is a tight distribution 
(standard deviation 0.3 dex), as also illustrated in panel (fe), where the PAH and FIR luminosities are plotted against each other. The 
meaning of the symbols is the same as in Figure |4] the solid line is the mean ratio assuming a slope of unity, and the dashed lines show the 
standard deviation from the mean. (The Hll-like ULIRG that lies well above the best linear PAH-FIR fit is F00397-1312; the PAHs in 
this system are corrected for significant extinction.) This ratio is very similar to that of the PG QSOs (Paper I), (d) L(PAH)/L(FIR) vs. 
effective silicate optical depth. The hne is a fit to HII/LINER ULIRGs. (F00397-1312, with the highest observed L(PAH)/L(FIR) value, 
is not visible in this plot). An anticorrelation exists between PAH/FIR ratio and silicate optical depth. 
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Fig. 21. — (a) Distribution of equivalent widt hs o f the 7.7 fim PAH feature for 1-Jy ULIRGs and average FIR-bright, FIR-faint, and FIR- 
undetected PG QSOs (see the legend in Figure [TH . Note the similarity between Seyfert 1 ULIRGs and QSOs, and the higher EW(PAH) 
values in HII galaxies and LINERs than in Seyferts. In the other panels, the 7.7 /im PAH equivalent width is plotted as a function of (b) 
the infrared luminosity, ( c) t he 30-to-15 fim flux ratio, and (d) the silicate optical depth. The meaning of the square (circle) symbols is 
the same as in Figure HI 1116 l l. These figures show that EW(PAH) is weakly correlated with infrared luminosity, strongly correlated with 
/30//15 (more so than with optical spectral type), and strongly correlated with silicate optical depth for HII/LINER, galaxies. The line in 
(d) is a fit to the HII/LINER detections. Seyferts deviate from the HII/LINER correlation and populate the lower left-hand corner of the 
plot. 
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Fig. 22.— (a) [Ne II] 12.8 fim luminosity, (6) [O IV] 25.9 ^lm luminosity, (c) [Ne V] 14.3 urn luminosity, and (d) [Ne VI] 7.65 fim 
luminosity vs. FIR luminosity of the 1-Jy ULIRGs and PG QSOs in the Spitzer sample as well as some optically-selected star bursts and 
Seyfert 2 galaxies observed with ISO. The meaning of the square and circle (other) symbols is the same as in Figure l4l I I12I I. The solid 
(dotted) line represents the mean (standard deviation) of the ratio L(emission line)/L(FIR) for HII ULIRGs. All ULIRGs lie close to the 
fitted [Ne II]-to-FIR ratio of IQ-^ ^^ but QSOs and Seyfert ULIRGs are increasingly above the values of the [O IV]- and [Ne V]-to-FIR 
luminosity ratios of HII ULIRGs. This result is readily explained if the AGN contributes increasingly to the fine-structure line emission at 
higher ionization levels. 
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Fig. 23. — [Ne III] 15.5 /im/[Ne II] 12.8 flux ratios vs. (a) total infrared luminosity, (6) FIR luminosity, (c) 25-to-60 flux ratio, (d) 
30-to-15 um flux rati o, and (e) 7.7 fim PAH equivalent width. The meaning of the small square and circle (other) symbols is the same as in 
Figure l4l I I12I and ll6l l. In addition, the small, medium-size, and large black circles in panel (e) correspond to the average FIR-undetected, 
FIR-faint, and FIR-bright PG QSOs from Paper II. No obvious trend is seen with FIR or total infrared luminosity among ULIRGs, but a 
clear dependence is seen on optical spectral type and IR continuum colors. 
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Fig. 24. — [O IV] 25.9 /im/[Ne II] 12.8 ^im flux ratios vs. (a) total infrared luminosity, (fe) FIR luminosity, (c) 25-to-60 /^m flux ratio, (d) 
30-to-15 um flux ratio, and (e) 7.7 /^m PAH equivalent width. The meaning of the small square and circle (other) symbols is the same as in 
Figure l4l Ill2l and ll6l l. This ratio among ULIRGs is slightly larger on average among ULIRGs with larger far- or total infrared luminosity. 
[O IV]/[Ne II] is more clearly larger among warmer, Seyfert ULIRGs. PG QSOs nicely fit along this excitation sequence. 
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Fig. 25. — [Ne V] 14.3 /^m/[Ne II] 12.8 flux ratios vs. (a) total infrared luminosity, (fe) FIR luminosity, (c) 25-to-60 flux ratio, 
(d) 30-to-15 mum fl ux r atio , an d (e) 7.7 fim PAH equivalent width. The meaning of the small square and circle (other) symbols is the 
same as in Figure l4l I ll2l and ll6l l. This ratio among ULIRGs is slightly larger on average among ULIRGs with larger FIR or total infrared 
luminosity. [Ne V]/[Ne II] is more clearly larger among warmer, Seyfert ULIRGs. PG QSOs nicely fit along this excitation sequence. 
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Fig. 26.— [Ne V] 14.3 fim/[Ne II] 12.8 fim vs. [O IV] 25.9 Atm/[Ne II] 12.8 fim for ULIRGs and PG QSOs in our sample. The meaning 
of the square and circle (other) symbols is the same as in Figure 141 l|12| l. A clear positive correlation is observed, representing an excitation 
sequence anchored with the HII-like/LINER ULIRGs, moving up in excitation level to the Seyfert 2 and Seyfert 1 ULIRGs, and ending 
with the PG QSOs. The solid diagonal line in these diagrams indicates constant [Ne V]/[0 IV] and changing [Ne II]. This line may be 
interpreted as a mixing line if [Ne V] and [O IV] are only produced by an AGN and [Ne II] by starburst activity. The tickmarks along 
the line indicate the percent contribution of the starburst to [Ne II] . The anchor point of this line has some systematic uncertainty due to 
variable AGN physical conditions. 
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Fig. 27. — [Ne VI] 7.65 /im/[Ne II] 12.8 /^m flux ratios vs. (a) total infrared luminosity, (6) FIR luminosity, (c) 25-to-60 flux ratio, 
(d) 30-to-15 fi m flu x ratio, and (e) 7.7 /xm PAH equivalent width. The meaning of the square and circle (triangle) symbols is the same 
as in Figure l4l II12I I. [Ne VI] is detected only in Seyfert ULIRGs and QSOs. The small number of detections prevent us from looking for 
statistically significant correlations with infrared or FIR luminosity, optical spectral type, or continuum colors. 
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Fig. 28.— (a) [Ne VI] 7.65 fim/[Ne II] 12.8 ^im vs. [Ne VI] 7.65 ^im/[0 IV] 25.9 ^im and (fc) [Ne VI] 7.65 fim/[Ne III] 15.5 fim vs. [Ne VI] 
7.65 ^m/[Ne V ] 14 .3 fim for ULIRGs and PG QSOs in our sample. The meaning of the square and circle (triangle) symbols is the same 
as in Figure l4l l|12| l. The grids show the predictions from the dusty, solar-metallicty, radiation pressure-dominated photoionization models 
of Groves et al. (2004) for two different power-law slopes [fi, oc u", a = —1.2 (blue) and a = —2.0 (red)]. The two parameters that vary 
across the grids are the ionization parameter (log U = —4 - 0) and the hydrogen density (log nj{ = 2-4). The color hues correspond to 
low values (light color) and high values (dark color); e.g., dark blue means high density or high ionization parameter. Poor agreement is 
seen in (a) due to starburst contamination of [Ne II] or problems with the model (highlighted by disagreements in both axes). Panel (b) is 
more useful, since it uses only high ionization lines of a single element (Ne). The good agreement means that AGN may power not only 
[Nc VI] and [Ne V], but also [Ne III], in these systems. 
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Fig. 29. — Luminosity of H2 S(l) 17.04 ^m vs. various star formation indicators: (a) FIR luminosity, (6) luminosity of [Ne II] 12.8 ^m, 
and (c) luminosity of PAH 7.7 ^m. The meaning of the square and circle (other) symbols is the same as in Figure 141 l|12|l . The solid and 
dotted lines show the mean and standard deviation of L[H2 S(1)]/L(X) (where X = FIR, etc.) for HII ULIRGs. A positive correlation is 
seen between these quantities for both ULIRGs and PG QSOs. However, Seyfert ULIRGs and especially PG QSOs have larger H2-to-FIR, 
H2-to-[Ne II], and H2-to-PAH ratios than HII ULIRGs, suggesting that H2 emission is influenced by the presence of the AGN in these 
objects. The LINERs also have significantly larger ratios and dispersions in panels (b) and (c). 
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Fig. 30. — Histograms of H2 (a) S(2)/S{1), (b) S(3)/S(l), and (c) S(3)/S(2) flux ratios (uncorrected for extinction) as a function of 
/25//6O in ULIRGs and QSOs. The ULIRG distributions are divided at the median value of /25//60- These ratios are indicative of the 
excitation temperature of H2 gas, but can also be affected by extinction and variable ortho-to-para ratios. Warm ULIRGs have smaller 
S(2)/S(l) ratios on average than cool ULIRGs by 0.1 dex, with a K-S (Kuipcr) probability of arising from the same distribution of <0.1% 
(1%). Warm ULIRGs have larger S{3)/S(2) ratios on average than cool ULIRGs by 0.15 dex, with a K-S (Kuiper) probability of arising 
from the same distribution of <0.1% (2%). Warm and cool ULIRGs are statistically indistinguishable in S(3)/S(l). 
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Fig. 31. — H2 excitation diagram for a few ULIRGs in our sample. The excitation temperature (Tex) is the reciprocal of the slope 
between any two data points. Solid lines join points where both lines are detected, and dashed lines join points where at least one line is 
an upper limit. Black lines join immediately adjacent lines, and red lines join the the S(l) and S(3) transitions. Error bars are 2a. 
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Fig. 32. — Difference between the excitation temperatures derived from the H2 S{3)/S(2) and S(2)/S(l) flux ratios vs. (a) /25//6O a-nd 
(b) the 9.7 /^m siHcate effective optical depth. The meaning of the symbols is the same as in Figure m The size of the symbol reflects the 
relative uncertainties on each data point, where the quartile of most certain points are the largest and the quartile of least certain points 
are the smallest. All objects should be above the solid line, unless extinction and/or ortho-para effects are at play. The small downward 
arrow on the right in each diagram reflects the effect of changing the ortho-to-para ratio from 3 to 2, while the long arrow reflects the 
effect of an extinction Ay = 10. Signiflcant trends among ULIRGs are seen, with decreasing temperature difference with increasing silicate 
optical depth and decreasing /25//60- ULIRGs with extinction greater than the median (and /25//6O lower than the median) have a lower 
temperature difference by 60 — 70 K, with a K-S (Kuipcr) signiflcance of <0.1% (6—7%). The trend with extinction implies extinction of 
the molecular lines. When corrected for this extinction, the trend with /25//6O will lessen or disappear. 
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Fig. 33. — (a) Average IRS spectrum of PAH-dominated ULIRGs from Figure lTOl minus one object with high obscuration (F00397— 1312). 
The average spectrum is from individual spectra normalized at 15 /im. The zoomed-in insert shows the detection of Hu« 12.4 fim. (b) 
[Ne III] 15.5 ^im/Hua and [Ne II] 12.8 /xm/Hua ratios derived from this average spectrum. The neon abundance relative to hydrogen 
increases from the lower-left portion of this diagram to the upper-right, as indicated by the solid iso-metallicity curves. The line ratios of 
the average spectrum suggest a neon abundance ~2.9 X solar, based on the Asplund et al. (2004) normalization (see text for a discussion 
of the uncertainties on this value) . 
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Fig. 34.— AGN/HII/PDR mixing diagram based on the Laurent et al. (2000) method, as modified by Armus et al. (2007): PAH (6.2 
/im) to continuum (5.3 - 5.8 fim) flux ratios vs. the continuum (14 - 16 fim) / (5.3 - 5.8 /im) flux ratios. The meaning of the symbols is 
the same as in Figure |4] The zero points for the pure HII region (upper-right) and PDR (lower-right) are from Armus et al. (2007) and the 
zero point for the pure AGN (lower-left) corresponds to the average value for the FIR-undetected PC QSOs to reduce possible starburst 
contributions to the continuum emission (Paper II). Note that the percentages included here are percentages of the 5.3 — 5.8 /im continuum; 
actual AGN fractional contributions to the bolomotric luminosities will be lower. 
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Fig. 35. — Comparison of AGN fractional contributions to the bolometric luminosities of ULIRGs and PC QSOs derived from 6 different 
methods: the {[O IV] 25.9 /xm/[Ne II] 12.8 ^lm) ratio, the ([Ne V] 14.3 ^m/[Ne II] 12.8 urn) ratio, the equivalent width of PAH 7.7 urn, 
the PAH (5.9 - 6.8 fim) to continuum (5.1 - 6.8 fim) flux ratio combined with the continuum (14 - 15 fim) / (5.1 - 5.8 /im) flux ratio, 
the MIR blackbody to FIR flux ratio, and the /30//15 continuum flux ratio. Squares are actual AGN fractional contributions, while small 
circles are upper limits on the AGN fractional contributions derived from the flne structure line ratios. The colors represent deviations 
(x — y) from the line of equalit y, on a logarithmic scale. Good agreement is seen on average between the various methods. See Appendix 
A, Section 7.1.1, and Tables [9l - I 121 for a description of each method and the results of these comparisons. 
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Fig. 36. — AGN fractional contributions to the bolometric luminosities of ULIRGs and PG QSOs, averaged over all 6 methods, vs. (a) 
optical spectral type; (b) infrared luminosity; (c) /25/ /60; (d) nuclear separation; and (e) interaction class. The meaning of the small 
square (large circle) symbols is the same as in Figure l4l III6I I. Overlaid crosses indicate ULIRGs with higher than average MIR extinction. 
The average (±1 standard error) AGN contributions in each horizontal bin are connected by the solid (dashed) lines. These binned points 
include ULIRGs only. The scatter in spectral type and interaction class [panels (a) and (e)] is added artificially for clarity. Strong positive 
correlations are detected between AGN contribution and optical spectral type / continuum slope: more Seyfert-like and warmer galaxies 
are increasingly AGN-dominated. Weaker, though significant, trends are seen for the other 3 independent variables: AGN contribution is 
highest at the highest infrared luminosities, smallest nuclear separations, and latest interaction classes. In all cases but infrared luminosity, 
the addition of PG QSOs cither extends or enhances these trends. 
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Fig. 37. — Effective optical depth of the 9.7 fim silicate feature vs. equivalent width of the 7.7 /im PAH feature as a function of (a) AGN 
fractional contribution to bolometric luminosity, (b) nuclear separation, and (c) interaction class. The meaning of the symbols is the same 
as in Figure |4] For clarity, we divide this diagram into three regions of roughly equal numbers. AGN contribution increases as one moves 
from Region 1 through Region 2 and into Region 3. There is also a trend toward smal ler nuclear separations and later merger stage along 
the sequence Rl — R2 — R3, but it is a weak trend with significant scatter (see Table [T3t . 
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Fig. 38. — Eddington ratio from photometry vs. (a) optical spectral type; (b) infrared luminosity; (c) /2 5// 60; (d) nuclear separation; 
and (e) interaction class. The meaning of the small square (large circle) symbols is the same as in Figure |4] CSJ . Overlaid crosses indicate 
ULIRGs with higher than average MIR extinction. The average (±1 standard error) Eddington ratios in each horizontal bin are connected 
by the solid (dashed) lines. These binned points include ULIRGs only. The scatter in spectral type and interaction class [panels (a) and 
(e)] is added artificially for clarity. There are no obvious correlations between Eddington ratios and any of these paramters. The Eddington 
ratios appear to be highest at the smallest nuclear separations and latest interaction classes, and the addition of PC QSOs seems to either 
extend or enhance these trends, but these results cannot be confirmed in a more rigorous statistical analysis. 
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Fig. 39. — Eddington ratio from photometry vs. (a) optical spectral type; (b) infrared luminosity; (c) /2 5// 60; (d) nuclear separation; 
and (e) interaction class. The meaning of the small square (large circle) symbols is the same as in Figure |4] CSJ . Overlaid crosses indicate 
ULIRGs with higher than average MIR extinction. The average (±1 standard error) Eddington ratios in each horizontal bin are connected 
by the solid (dashed) lines. These binned points include ULIRGs only. The scatter in spectral type and interaction class [panels (a) and 
(e)] is added artificially for clarity. There are no obvious correlations between Eddington ratios and any of these paramters. The Eddington 
ratios appear to be highest at the smallest nuclear separations and latest interaction classes, and the addition of PC QSOs seems to either 
extend or enhance these trends, but these results cannot be confirmed in a more rigorous statistical analysis. 



